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ABSTRACT 

We present an analysis of the 2 — 10 keV X-ray emission associated with the active 
galactic nuclei (AGNs) in brightest cluster galaxies (BCGs). Our sample consists of 32 
BCGs that lie in highly X-ray luminous cluster of galaxies (Lx-ray (0.1 — 2.4 keV) > 
3 X 10"*^ erg s~^ ) in which AGN-jettcd outflows are creating and sustaining clear X- 
ray cavities. Our sample covers the rcdshift range < z < 0.6 and reveals strong 
evolution in the nuclear X-ray luminosities, such that the black holes in these systems 
have become on average at least 10 times fainter over the last 5 Gyrs. Mindful of 
the potential selection effects that may affect our results, we propose two possible 
scenarios to explain our results: 1) either that the AGNs in BCGs with X-ray cavities 
are steadily becoming fainter, or more likely, 2) that the fraction of these BCGs with 
radiatively efficient nuclei is decreasing with time from roughly 60 per cent at z w 0.6 
to 30 per cent at z « 0.1. Based on this strong evolution, we predict that a significant 
fraction of BCGs in z « 1 clusters may host quasars at their centres, potentially 
complicating the search for such clusters at high redshift. In analogy with black-hole 
binaries and based on the observed Eddington ratios of our sources, we further propose 
that the evolving AGN population in BCGs with X-ray cavities may be transiting from 
a canonical low/hard state analogous to that of X-ray binaries to a quiescent state 
over the last 5 Gyrs. 

Key words: Galaxies: clusters: general - X-rays: galaxies: clusters - cooling flows - 
galaxies: jets - black-hole physics - accretion, accretion discs 



1 INTRODUCTION 

The interplay between the accretion of material onto a su- 
permassive black hole (SMBH) and the release of energy 
through radiation and outflows is known as active galac- 
tic nucleus (AGN) feedback. Some of the strongest cases of 
AGN feedback are seen in clusters of galaxies, where the cen- 
tral AGN hosted by the brightest cluster galaxy (BCG) is 
capable of driving large jetted outflows flUed with relativis- 
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tic plasma. As the jetted outflows propagate through the 
intracluster plasma, they push aside the hot X-ray emitting 
gas, creating cavities that are detectable as regions of re- 
duced surface brightness in X-ray images. These X-ray cav- 
ities act as calorimeters and provide a unique opportunity 
to directly measure the work done by the AGN on the sur- 
rounding medium. 

X-ray cavities are therefore extremely useful tools for 
studying the details of AGN feedback. Studies at low red- 
shifts have shown that these structures are predominantly 
found in c ool-core clusters of ga, l axies (detectio n rate > 90 
per cent; iDunn fc Fabiaiil l2006l : iFabianI |2012|) . Cool-core 
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clusters have highly peaked X-ray surface brightness pro- 
files, and short central cooling times that are often shorter 
than the Hubble time. The hot X-ray gas at the centres 
of these clusters should therefore have had the time to 
cool, and large flows of cooling material, known as cooling 
flows, would naiv ely be expected in the central regions of 
these objects fsee ISarazinlll98^ : iFabian et al.lll984l : FabianI 
1 1994 ). However, both Chandra and XMM observations 
have shown that there is significantly less cooling material 
at the centres of these clusters than expected from standard 
cooling-flow models llBohringer et aL 200 ll: Tamura et al.l 



20011: jPeterson et al.l l200ll l2003l : iPeterson fc FabianI l2006l : 

03: 



McNamara fc NulsenI 20071 ). This is known as the cooling- 



flow problem, and feedback from the central AGN is thought 
to be the leading mechanism that prevents the hot X-ray gas 
from cooling, by inflating the X-ray cavities and propagating 
energy through shock and sounds waves. Studies of nearby 
clusters {z < 0.3) have shown that the energetics of these 
X-ray cavities are indeed capable on average of preventin; 
the hot X-ray gas from coohng j Birzan et all |2004 |200 
Dunn et al.ll2005l: JPunn fc FabianI 
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2007 ; iDunn ct al."2010'; 'C a.vagnolo et al.ll201ol : jPong et Zl 
201C ; lo^Sullivan ct al.. 
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Recently, we have extended the sample of known X- 
ray cavities into the higher-red shift Universe (0.3 < z < 0.6; 
iHlavacek-Larrondo et aLll2012l . h ereafter HL201 2) usii ig the 
MAs s ive Cluster Survey (MA CS; lEbeling et al]|200ll . |2007| . 
I2OIOI : iMann fc Ebeling||2012l ). The MACS survey compiled 
the flrst large sample of very X-ray luminous clusters of 
galaxies at intermediate to high redshift, and consists of 124 
spectroscopically confirmed clusters at 0.3 ^ 2 ^ 0.7 (see 
Fig. [ij . Our work on the MACS clusters showed that X-ray 
cavities remain common in high-redshift cool-core clusters 
and that the energetics are still capable of preventing the 
surrounding gas from cooling. In particular, we found no ev- 
idence for evolution in any of the jetted outflow properties, 
implying that extreme mechanical AGN feedback has been 
in place for at least the past 5 Gyrs. We also noted that many 
of the clusters with AGN outflows in the MACS sample had 
bright X-ray AGN. This is ra rely seen in clusters of similar 
lumin osities at low redshift (|Hlavacek-Larrondo fc FabianI 
I2OIII . hereafter HL2011), where most of the X-ray lumi- 
nous clusters ai z < 0.3 have no detectable X-ray nucleus, 
suggesting that we may be seeing some form of evolution in 
the radiative properties of the black holes in these BGGs. 
AGN feedback appears to be operating differently in high- 
redshift BGGs with X-ray cavities, at least in terms of the 
radiative properties of the AGN. 

Such an evolution is indicative of the predicted evo- 
lution between "quasar-mode" feedback and "radio-mode" 
feedback. Local active SMBHs are mostly thought to be in 
radio-mode feedback, where the accretion rates are low and 
the black hole is capable of driving powerful jetted outflows. 
This explains the behaviour of nearby galaxy cluster cores 
{z^O.2 — 0.3), which appear to have radiatively inefficient 
central AGN, yet are capable of driving large kpc scale jet- 
ted outflows. On the other hand, quasar-mode feedback is 
powered by a radiatively efficient black hole accreting at 
rates near the Eddington limit. This feedback mode has been 
invoked to explain the apparent deficit of extremely lumi- 
nous galaxies in the galaxy luminosity function as well as 
the tight relations between black hole mass and host galaxy 
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Figure 1. Cluster X-ray luminosity versus redsh ift. Shown are 
the extend ed Brightest Cluster Survey (eBCS: lEbeling et al.l 
Il998l.l2000l ). the 400 deg^ cluster s urvey jBurenin et ahlpOOTl ). the 
XMM — newton cluster survey jMehrtens et al.l 12013) " and the 
MACS survey which finds the luminous clusters at 0.3 < z < 0.7. 
The black diamonds highlight the clusters in our sample and the 
horizontal dashed line shows our > 3 X 10** erg s~^ luminosity 
cut. The empty black diamonds show the clusters not part of 
eBCS. Here, we use the cluster luminosities from Table|2]to illus- 
trate them on the plot. Note that, we do not include 4C-I-55.16 
and H182H~643 since they are not part of any catalogue. 



properties (e.g. stellar bulge mass and stellar velocity dis- 
persion). Simulations have also demonstrated that quasar- 
mode feedback can regulate star formation, by heating 



star formation (iDi Matteo et al.l 


2OO5I: iCroton ct al.' '20061: 


Bower et al. 2006: Soringel et al. 


I2OO5I: ISiiacki & Spring's! 


20061: iMerloni fc Heinzll2008l). Although the quasar luminos- 



ity fu nction peaks at high redshift (z « 2 — 3, e.g. lWall et al.l 
l2005h . indicating that AGN in the past are more radiatively 
efficient, and that present day SMBHs are predominantly 
operating in a radiatively inefficient mode, it is still not clear 
how AGN transit from one mode to the other. Clusters of 
galaxies provide one of the most direct pieces of evidence 
for AGN feedback, with clear evidence that black holes can 
have a substantial impact on their surrounding medium (i.e. 
by the observations of X-ray cavities excavated by the AGN- 
driven outflow). By studying the radiative properties of the 
central AGN in clusters with known X-ray cavities (i.e. with 
known radio-mode feedback taking place) as a function of 
redshift, we can therefore provide insight into how quasar- 
mode and radio-mode feedback evolve with cosmic time. 

The objective of this paper is therefore to further anal- 
yse the radiative properties of BGGs with X-ray cavities. 
We flrst describe the details of our sample selection in Sec- 
tion 2, and derive relevant cluster properties in Section 3. 
In Section 4, we measure the nuclear X-ray luminosities, 
and then radio luminosities in Section 5. Section 6 presents 
the results and Section 7 discusses high-redshift clusters. 
Finally, in Section 8 we discuss the implications of our re- 
sults and in Section 9, we present the conclusions. We adopt 
Ho = 70 km s"^ Mpc"^ with = 0.3, = 0.7 through- 
out this paper. All errors are la unless otherwise noted. 
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2 SAMPLE SELECTION 

Our aim is to study the radiative evolution of the AGN in 
BCGs with known X-ray cavities. We begin by selecting the 
intermediate to high redshift clusters (0.3 < z < 0.6) with 
known X-ray cavities, and base this on our previous work in 
HL2012, where we identified X-ray cavities in MACS using 
Chandra images. Chandra observations are required since 
it is the only X-ray telescope that has the resolving power to 
identify X-ray cavities beyond the local Universe {z^Q.Q2). 

76 MACS clusters have archival Chandra observa- 
tions. For each of these, we computed unsharp-masked and 
elliptical-subtracted X-ray images. Both techniques are used 
to enhance deviations in the original images. Using these 
processed images, we searched for X-ray cavities, and found 
13 clusters with clear cavities (which could be clearly seen 
in the original images and processed images), as well as 7 
with potential cavities (which could only be clearly seen in 
the processes images). Deeper observations are needed to 
confirm the potential cavities. We therefore only consider 
the first 13 clusters as our high-redshift {z > 0.3) sample of 
objects with clear evidence for active mechanical feedback 
from their central AGN. These are listed in Table [T] 

In order to study the radiative evolution of AGN feed- 
back with cosmic time, we need to identify a comparison 
sample of low-redshift clusters with similar properties to 
the MACS clusters. We begin by searching the literature 
for clusters with known X-ray cavities at z < 0.3. In total, 
we find 40 clusters with reported X-ray cavities. 

The MACS survey is based on an X-ray fiux-limited 
sample, and only comprises highly X-ray luminous clus- 
ters. The fiux limit for MACS is 10"^^ erg cm"^ s"^ in the 
0.1 — 2.4 keV energy band. For our cosmology, this corre- 
sponds to a luminosity of 3 x 10*^ erg s^'^ at z — 0.3 (min- 
imum redshift for MACS). MACS therefore only consists of 
clusters that have luminosities larger than 3 x 10**^ erg s^^ . 
In order to compare this population of clusters to the low- 
redshift counterparts, we must only consider low-redshift 

clusters with luminosities above 3 x 10'*^ erg s~^ . This is 

1 \ 



Table 1. Sample of luminous clusters with known X-ray cavities. 



because X-ray luminous clusters ( Lx > 10''*"*^ erg s 
with strong cool cores, such as the MACS clusters with 
X-ray cavities, require extreme feedback from their cen- 
tral AGN to offset cooling of the surrounding medium 
> 10*** ^^ergs *), whereas less X-ray luminous 
clusters (« lO*'^ erg s~^ ), even with a strong cool core, only 
require some 10*'^ erg s^* of feedback from their central 
AGN. This translates to smaller and less powerful X-ray 
cavities in the less luminous clusters. If we were to include 
these in our sample, we would not be probing the same out- 
fiows as in MACS. Out of the 40 clusters with clear X-ray 
cavities at z < 0.3, we therefore only consider those that 
have similar luminosities to MACS ( Lx > 3 x 10** erg s~* 
in the 0.1 — 2.4 keV energy band). We use the available total 
cluster fiuxes in vari ous surveys based on the ROSAT AU- 
Skv Survev (fl^gg: lEbeling et all 1 19961 . Il998l . |2000| . |2002| : 



iBohringer et al. I I2OO4I ) to determine if a cluster is luminous 
enough to be included in our sample. If a source was present 
in more than one catalogue, we only considered the most 
recent catalogue. Our final sample is shown in Table [1] and 
examples in ascending order of redshift are shown in Fig. (2] 
The catalogue cluster fiuxes are also shown in Table [21 

Our sample also includes the powerful radio galaxy. 
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Cygnuf 


A 




0.0561 


Abell 1795 




U.UOo 


Abell 2029 




0.077 


Abell 2597 




0.085 


Abell 0478 




0.0881 


RXC J1558.3-1410 




u.uy / U 


RXC J1524.2-3154 




0.1028 


PKS 0745-19 




n 1 noQ 
U. iUZo 


Abell 2204 




0.1522 


Hercules A 




0.154 


Abell 0115 




0.1971 


ZwCl 2701 




0.2151 


MS0735.6-I-7421 


ZjWL^I io/U 


0.2160 


4C-I-55.16 






Abell 1835 




0.2532 


ZwCl 3146 




U.iyUb 


H182H-643 




0.299 


MACS 


J2140.2-2339 


Mb Zl6i -o-Zooo 


0.313 


MACS 


J0242. 5-2132 




0.314 


MACS 


J0547.0-3904 




0.319 


MACS 


J1931. 8-2634 




0.352 


MACS 


J0947. 2-1-7623 


RBS 0797 


0.354 


MACS 


J1532.8-(-3021 


RX J1532. 9-1-3021 


0.3613 


MACS 


J1720.2-I-3536 


Z8201 


0.3913 


MACS 


J0429. 6-0253 




0.397 


MACS 


J0159. 8-0849 




0.404 


MACS 


J2046. 0-3430 




0.423 


MACS 


J0913. 7+4056 


IRAS 09104-1-4109 


0.442 


MACS 


J1411. 3-1-5212 


3C295 


0.460 


MACS 


J1423. 8-1- 2404 




0.5449 



4C+55.16 (2 = 0.2412). 4C-I-55.16 is not part of any of the 
catalogues, but is embedded in a luminous cluster that has 
an X-ray luminosity within 200 kpc larger than many of the 
other clusters in our sample (see Table [JJ . It is therefore 
reasonable to assume that if it were part of the catalogues, 
its X-ray luminosity would be larger than 3 x 10** erg s~* , 
and therefore should be included in our sample. 

Finally, we include H1821-I-643 {z = 0.299), which 
is a quasar embed ded in a luminous cool-core cluster 
(|RusseU et aLlbOlOal V Since H1821+643 shows evidence for 
X-ray cavities that coincide with radio lobes, we included 
this object in the initial sample of 40 clusters with known 
X-ray cavities. H1821-I-643 is however not part of any of the 
catalogue s, and this quasar is heavily affected by pile-up. 
However. iRussell et al.l (|2010al ) have accounted for the pile- 
up by simulating the quasar PSF to disentangle the quasar 
emission from the intrac luster medium. We t herefore use 
the quantities quoted in I Russell et al.l (|2010a ) throughout 
this pa per instead of reanalyzing the source. Russell et al.l 
(|2010ah furthermore estimate that the cluster luminosity in 
the 0.5 — 7 keV range is ~ 10*^ erg s^* , of which most orig- 



inates from the 0.1 — 2.4 keV energy band. It is therefor 
reasonable to assume that this cluster is luminous enough 
to be included in our sample. 
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Figure 2. Examples of 4 clusters of galaxies in our sample in ascending order of redshift from left to right. The top panels show the 
0.5 — 7 keV images, whereas the bottom panels show the 3 — 7 keV images where the non-thermal nuclear emission from the central 
point source should start dominating the emission if present. As seen, the emission from the central point source starts to dominate over 
the extended cluster emission with increasing redshift in the bottom panels. For the 2 low-redshift clusters (Abell 0085 and Hercules A), 
most of the emission seen is extended thermal emission. 



3 CLUSTER X-RAY LUMINOSITIES 

Our final sample comprises 32 clusters of galaxies in the 
0.0 < z < 0.6 redshift range. It includes some of the well 
known lumino us clusters such as the Perseus cluster, MACS 
J1931. 8-2634 (|Ehlert et al.ll201lh and MS 0735. 6H- 7421, the 
latter of which has the most po werful mechanical out- 
burst known to date (10""^ erg s"^ ; iMcNamara et al.ll2005l : 
iGitti et al.|[2007l ). 

The cluster X-ray luminosities as derived from each of 
the catalogues are shown in column 4 of Table [2] and the 
catalogue references in column 8. As mentioned in the pre- 
vious section, H1821 -I-643 and 4C- I- 55.16 a re not part of the 
catalogues. However. lRussell et al.l (|2010al ) estimate that the 
cluster luminosity for H1821-I-643 in the 0.5 — 7 keV range 
is roughly 10*^ erg s~^ . According to their spectra, most of 
this emission originates from the 0.1 — 2.4 keV energy band, 
and we therefore approximate that the 0.1 — 2.4 keV cluster 
luminosity for this source is lO"'^ erg s"^ . For 4C+55.16, we 
compute a rough estimate of the cluster luminosity by fitting 
MEKAL model to the 0.5 — 7 keV Chandra emission within 
the ce ntral 200 kpc. The Ga lactic absorption is kept frozen 
at the lKalberla et all (|2005l ') value and the background was 
chosen to be located on the same detector but far from any 
cluster emission. 

We now focus on deriving the cooling luminosities in or- 
der to measure the power of the predicted cooling flows. We 
deflne the cooling luminosity as the 0.1 — 2.4 keV luminosity 
within the cooling radius, an d we define the cooli ng radius 
using the same definition as iRaffertv et all (|2006l '). i.e. the 
radius at which the cooling time (icooi) is equal to the 2 = 1 
look-back time. We adopt the same definition so that we can 
directly compare our results with theirs. For our cosmology, 
this corresponds to t^oox = 7.7 Gyrs. To determine the ra- 
dius at which tcooi = 7.7 Gyrs, we have computed detailed 
cooling time profiles for all of our clusters. We use the same 
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10 100 
Radius (kpc) 

Figure 3. Cooling time profiles for all of our clusters shown in 
Table [l] The c ooling time profile for H1821-I-643 was taken from 
I Russell et al.l ll2010al ). Deprojected thermal plasma parameters 
were used to derive the profiles. Each colour illustrates a different 
cluster. 

method as in HL2012, where fcooi is calculated with Eq. [l] 
Here rio is the electron density, kT is the gas temperature, 
Lx is the gas X-ray luminosity and V is the gas volume 
contained within a shell. 

5 1.92 no fcT y ... 
icool = z (IJ 

For each cluster, we selected the deepest observational 
Chandra data set available. These data were then processed. 
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Table 2. Cluster properties - (1) Name; (2) Chandra observation identification number; (3) Remaining exposure time of the Chandra 
observations after data reduction; (4) X-ray luminosity in the 0.1 — 2.4 keV energy band as derived from catalogues, see Column 8; 
(5) Cooling radius, defined as the radius where the cooling time equals 7.7 Gyrs; (6) X-ray luminosity within the cooling radius in 
the 0.1 — 2.4 keV energy band; (7) Cavity enthalpy (PVcavitios) taken from the lite rature; (8) Reference f or th e catalogue where th e 
X- ray luminosity in Co lum n 4 was taken from: i) E beling et al. in prepa ration.; ii) iBohringer et ah 1 liool); iii) lEbeling et al. 1 1 I2OO2I) : 
iv) lEbeling et al ] l|2000l 'l: v'l lEbeling et al 1 lliggsll: yil lEbeling et al.l " Cluster luminosity derived as the 0.1 — 2.4 keV luminosity 

within 200 kpc. ''Values derived from iRusselTetaL I ll2010al) . 
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(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


Name 


ObsID 


Exp. 


-^X , catalogue 


'^cool 


-^cool 


^^avities 


Reference 






(ks) 


(lO'^" erg s-l ) 


(kpc) 


(10'*'' erg s-i ) 


(10^58 erg) 




Perseus cluster 


4952 


149.9 


7.03 


110 


on A 1 c\ r\r\'~i 

3.394+0.003 


5.3 


(ill) 


Abell 0085 


904 


37.1 


5.2 


68 


1.00+0.01 


1.2 


(vi) 


Cygnus A 


1707 


9.2 


4.0 


59 


1.00+0.02 


84 


(ill) 


Abell 1795 


493-1-494 


36.3 


6.5 


89 


2.40+0.01 


4.7 


(v) 


Abell 2029 


891 


19.9 


8.9 


81 


3.17+0.01 


4.8 


(v) 


A 1 11 c\'~T 

Abell 2597 


7329 


59.1 


4.7 


92 


c\ on 1 n n 1 

2.39+0.01 


3.6 


(vi) 


Al 11 r\ A f^iTi 

Abell 0478 


1669 


40.3 


7.7 


95 


K .4 1 n nn 

5.34+0.02 


1.5 


(v) 


KXO J1558.J— 1410 


9402 


36.4 


3.4 


81 


1.46+0.01 




(11) 


RXC J1524.2— 3154 


9401 


41.5 


3.0 


78 


1 nn 1 n n 1 

1.89+0.01 




(11) 


PKb 0745-19 


12881 


117.6 


12.3 


101 


^ ^ ^ In nn 

7.74+0.02 


69 


(iii) 


Al 11 c\c\r\ A 

Abell 2204 


7940 


73.1 


13.3 


90 


^7 nr? 1 n no 

7.97+0.03 


4.0 


(vi) 


Hercules A 


5796-1-6257 


96.2 


3.6 


60 


n 1 n n 1 

0.63+0.01 


31 


(v) 


Ai_ 11 rMir" 

Abell 0115 


3233 


43.1 


10.0 


80 


1 00 1 n nn 

1.38+0.02 


43 


(v) 


ZwOl 2701 


3195 


25.9 


7.6 


79 


-i m 1 n no 

1.91+0.03 


350 


(v) 


Mb 0735.6+7421 


10470 


135.0 


5.46 


79 


-i nn 1 n nn 

1.99+0.02 


1600 


(iv) 


40-I-55.16 


4940 


73.6 


w 4.8° 


91 


nn 1 n no 

3.09+0.03 


10 




Ai_ 11 loor" 

Abell 1835 


6880 


113.8 


28.7 


100 


10.8+0.1 


47 


(v) 


ZwCl 3146 


909-F9371 


77.8 


20.6 


131 


12.8+0.1 


380 


(y) 


H1821-I-643 


9398-h9845+9846-F9848 


85 


w lO** 


w 90*' 


?» 10** 


14 




MACS J2140.2-2339 


4974-h5250+928 


76.8 


7.7 


107 


8.13+0.05 


3.4 


(i) 


MACS J0242. 5-2132 


3266 


8.6 


10.1 


110 


9.51+0.2 


9.0 


(i) 


MACS J0547.0-3904 


3273 


19.2 


5.9 


100 


3.1+0.1 


3.2 


(i) 


MACS J1931. 8-2634 


9382 


95.0 


12.6 


112 


3.47+0.03 


83 


(i) 


MACS J0947.2-F7623 


7902 


38.8 


13.7 


120 


15.3+0.1 


145 


(i) 


MACS J1532. 8-1-3021 


1649 


9.5 


11.8 


115 


8.0+0.1 


32 


(i) 


MACS J1720.2-I-3536 


3280+6107-^7718 


51.7 


9.7 


100 


4.9+0.1 


18 


(i) 


MACS J0429.6-0253 


3271 


21.4 


9.5 


105 


6.5+0.2 


2.2 


(i) 


MACS J0159.8-0849 


3265+6106+9376 


61.4 


11.5 


86 


6.0+0.1 


13 


(i) 


MACS J2046.0-3430 


9377 


35.9 


6.9 


81 


5.0+0.1 


23 


(i) 


MACS J0913.7-h4056 


10445 


70.4 


9.7 


107 


8.8+0.1 


150 


(i) 


MACS J1411. 3-1-5212 


2254 


76.4 


6.7 


60 


2.4+0.1 


49 


(i) 


MACS J1423.8-I-2404 


4195 


106.9 


9.4 


95 


9.2+0.1 


61 


(i) 



cleaned and calibrated using the latest version of the CIAO 
software (ciAOv4.4, CALDB4.4.9/10), and starting from the 
level 1 event file. We applied both CTI (charge time inter- 
val) and time-dependent gain corrections, as well as removed 
flares using a 3cr threshold. We also exclude point sources 
during the fit, since these are more likely dominated by 
non-thermal emission. When a cluster was observed multiple 
times with the target centred on the same detector and in 
the same observing mode (faint or vfaint), we combined 
the different observations only if this would improve the im- 
age quality significantly. However, for MS 0735.6+7421 and 
the Perseus cluster, we only consider the deepest observa- 
tion since they are sufficient to obtain detailed cooling time 
profiles. 

The thermal gas properties are then determined by se- 
lecting a set of annuli containing the same number of counts 
and centred on the X-ray peak. The total number of annuli 



depend on the data quality, but we chose the minimum num- 
ber of counts so that we have at least 5 data points in the 
cooling time profile within 100 kpc, allowing us to estimate 
accurately the cooling radius when fcooi = 7.7 Gyrs. For the 
more nearby and well studied objects with deep Chandra 
observations, we allowed several tens of thousand counts per 
annulus, but for the most distant and poorly observed ones, 
we could only allow some 2000 counts per annulus. The back- 
ground region was chosen either within the same chip but far 
from any cluster emission or, for the more nearby objects, 
within a chip located on the same detector. Since we focus 
on the central 200 kpc of our clusters, where the emission 
remains strong, the choice of background did not signifi- 
cantly affect our results. We then proceeded in deprojecting 
the data using the Direc t Spectral Deprojection method of 
ISanders fc FabianI (|2007l ). For each deprojected spectrum, 
we fitted an absorbed (Galactic) mekal model. In some an- 
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nuli, it was difficult to constrain the abundance. In this case, 
we kept the abundance value frozen at the cluster average 
which we determined by selecting a region within the en- 
tire cluster (r < 200 kpc). Point sources were excluded from 
the fitted region. The derived deprojected spectral quanti- 
ties were then used as an estimate of the plasma parameters, 
allowing us to derive tcooi based on Eq. [T] 

Once we obtained the cooling time profiles (Fig. [3]), we 
determined the cooling radius, defined as the radius when 
tcooi = 7.7 Gyrs. The cooling luminosity (Lcooi) then corre- 
sponds to the 0.1 — 2.4 keV luminosity within this radius. We 
show our results in Table [21 Column 6. For H1821+643, we 
use the cooling time profile derived bv lRussell et al.l (|2010al l 
to estimate the cooling radius when t^ooi = 7.7 Gyrs (~ 90 
kpc), and then estimate the cooling luminosity based on the 
av erage heating required within the cooling radius as derived 
bv lRussell et ai.l (|2010al . ; 10*^ erg s"^ kpc"^). 



4 NUCLEAR X-RAY LUMINOSITIES 

Chandra X-ray observations are used to isolate the nuclear 
X-ray emission for our sample of BCGs. Some of our BCGs 
are dominated by non thermal emission at X-ray wave- 
lengths, while others have no detectable X-ray nucleus (e.g. 
A1835, A478, Z2701). We therefore derive the nuclear lu- 
minosities using 2 methods. The first is used to estimate 
the fiuxes in a systematic way for all objects, regardless of 
the nuclear spectrum and is based only on the core count 
rate. The second, focusses on the objects with clear power- 
law components in their spectrum, and we determine their 
nuclear luminosities by fitting a non-thermal model to the 
spectra. The two methods are described in the following sec- 
tions. 



4.1 Nuclear luminosities based on a core count 
rate 

For the first method, we use the web interface of pimmsQ 
(jMukail 1 19931 ). which allows us to convert a count rate 
into a flux. We first consider the highest redshift object 
in our sample, MACS J1423.8-f 2404 located &t z = 0.5449 
(1" = 6.4kpc). We select a 2 by 2 pixels square region in 
the 0.5 — 7 keV Chandra image, which corresponds to a 
1" X 1" (or 6.4kpcx6.4kpc) square region. Then, we select 
a surrounding region as a background, and more precisely, 
we consider a square annulus centred on the same position, 
with an outer 6 by 6 pixels square and an inner 4 by 4 pixels 
square. The total count number of the core is then calcu- 
lated by subtracting the background emission scaled to the 
same number of pixels as the central 2 by 2 pixels square. 
PIMMS is then used to convert the count rate into a flux, and 
then a luminosity in the 2 — 10 keV energy band. A power- 
law with photon index of 1.9 is used as a model, but our 
results are not sensitive to our value of the index, at least 
within ±0.2. We then proceed with the other objects, and 
select a central 6.4kpcx6.4kpc square for each cluster. The 
background is selected in a similar way, as a surrounding 
square annulus. Since X-ray counts are governed by Poisson 

^ http://hcasarc.gsfc.nasa.gov/Tools/w3pimms.html 
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Figure 4. Logarithm of the 2 — 10 kcV X-ray core luminosity in 
erg s~^ as a function of redshift, derived from the observed count 
rates (i.e. not corrected for any internal absorption, if present). 
The Perseus cluster, Hercules A, H1821+643, Cygnus A and 3C 
295 are highlighted with the red stars. Note that the values for 
Perseus and H1821-I-643 have been taken from the literature and 
have therefore already been corrected for internal absorption. The 
blue points illustrate the sources with evidence of central non- 
thermal point source associated with the BCG. 



noise, we can estimate that the error associated with each 
square region is given as the square root of the number of 
counts. Using the error propagation equation, we then esti- 
mate the la noise level for the background-subtracted count 
rate which we show in Table [3] 

We consider the same region for all objects in terms of 
kpc'^, and not in terms of arcsec'^. This is because a 1" x 1" 
square region at high redshift is sampling more of the cluster 
emission than a 1" x 1" square region at low redshift. For 
the nearest objects, the central 6.4kpcx6.4kpc square region 
therefore extends well over a dozen pixels. We choose this 
particular dimension because we are limited by the Point 
Spread Function (PSF) of the telescope for the high-redshift 
objects. Chandra's PSF is on the order of r = 1", and within 
this radius, around 90 per cent of the energy of a point 
source would fall in. For a 1" x 1" square region, some 70 
per cent of the energy of a point source falls in. If we were 
to choose an even smaller region, we would be missing most 
of the nuclear fluxes for our high-redshift objects. Although 
this method does have its limitations, it should only affect 
the high-redshift objects by slightly underestimating their 
fluxes, and the low-redshift objects by slightly overestimated 
their fluxes since it will be including the counts from the 
nucleus and surrounding thermal emission. This means that 
any evolution we find, is probably a lower limit to the true 
evolution, which should be even more rapid. Our results are 
illustrated in Fig. |3] 

Note that, we could not derive a core flux for Perseus or 
H1821-I-643 based on the Chandra observations which are 
affected by pileup. Pileup occurs when two or more photons 
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Table 3. Core X-ray count rates (0.5 — 7 keV), as well as observed fluxes and luminosities in the 2 — 10 keV energy band within a 
6.4kpcx6.4kpe square region. Internal nuclear absorption at the redshift of the source has not been co nsidered here (see Table 4 for 
corrected l uminosities). "The X- ray nuclear luminosity for Perseus was taken from lMerloni fc Heind l|2007l '). ''Values for H1821-I-643 were 
taken from lRussell et aL I ll2010al') . 



(1) 


(2) 


(3) 


(4) 


Name 


Count rate 


Fx , core 


Lx.core 




(Cts s-i) 


(10-1* erg cm-2 s"! ) 


(10*2 erg s-i ) 


Perseus cluster"" 






25-17 


Abell 0085 


0.007+0.002 


4.3itl.2 


0.31±0.09 


Cygnus A 


0.307+0.007 


131±3 


9.8±0.2 


Abell 1795 


0.015+0.002 


6.5±0.8 


0.62±0.08 


Abell 2029 


0.032+0.003 


14.1±1.3 


2.1±0.2 


Abell 2597 


0.008+0.001 


3.4±0.5 


0.61±0.08 


A K/iU r\A'7Q 


U.UUOztU.UUi 


o.yitU.y 


n Q-Ln o 
U.oitU.Z 


rCAO JiOOo.o— i4iU 


U.UU00±U.UUUo 


2.3it0.3 


0.00±0. U8 


rtAO Jioz4.z— olo4 


n r\r\A -i-Ci nm 
U.UU4±U.UU1 


1.7ib0.5 


0.4it0.1 


FKb U745-iy 


U.OUOo+U.UUUo 


2.4±1.6 


0.6it0.4 


Abell 2zU4 


U.UU44+U.UUU/ 


3.1it0.5 


1.9ib0.3 


Hercules A 


U.UUzy+U.UUUvS 


1.4it0.2 


n m _l_n on 
0.yi±0.U9 


Abell Uiio 


U.UUU4+U.UUU2 


0.3it0.2 


0.3ib0.2 


ZjWLjI Z (Kjl 


n nnnc;_Ln nnnQ 
U . UUUoitU . UUUo 




n Q-Ln o 
(i.oztV.Z 


Mb U/v5o.bH- ( 4il 


U.UUUo±U.UUUl 




0.3ib0.1 


4C+oo.io 


U.UlUo ±U.UUUo 


4.9lt0.2 


8.6ib0.4 


Ai-.^ii lOQc; 
Abell loou 


U.UUUo+U.UUUo 


0.4it0.2 


0.8ib0.4 


r7,^,(^l n-\ Ad 

ZjWL/1 oi4b 


U.UU21 ±U.UUUo 


1.5it0.2 


3.9it0.6 


rlli521H-o4o 






4200_io(, 


MACS J2140. 2-2339 


0.0029+0.0003 


1.3±0.1 


4.3±0.5 


MACS J0242. 5-2132 


0.0026+0.0009 


1.8±0.6 


5.8±2.0 


MACS J0547.0-3904 


0.0037+0.0005 


2.5±0.3 


8.4±1.1 


MACS J1931. 8-2634 


0.0050+0.0002 


3.7±0.2 


15.3±0.7 


MACS J0947.2+7623 


0.0127+0.0006 


5.7±0.2 


24.0±1.1 


MACS J1532. 8+3021 


0.0003+0.0002 


0.15±0.09 


0.7±0.4 


MACS J1720. 2+3536 


0.0004+0.0002 


0.3±0.1 


1.6±0.6 


MACS J0429. 6-0253 


0.0011+0.0004 


0.8±0.3 


4.4±1.6 


MACS J0159. 8-0849 


0.0016+0.0003 


1.1±0.2 


6.2±1.1 


MACS J2046. 0-3430 


0.0007+0.0003 


0.5±0.2 


3.3±1.2 


MACS J0913. 7+4056 


0.0070+0.0004 


4.7±0.3 


33.6±2.0 


MACS J1411. 3+5212 


0.0038+0.0003 


2.5±0.2 


20.0±1.5 


MACS J1423. 8+2404 


0.0008+0.0002 


0.33±0.07 


3.9±0.8 



are detected as one event (see for more details lDavisll200ll : 
iRussell et alll2010bh . and the amount of pileup can be esti- 
mated by comparing the fraction of good (grades 0,2,3,4,6) 
to bad grades (grades 1,5,7) for each point source. Typi- 
cally, pileup becomes problematic when the fraction of bad 
grades exceeds 10 per cent of the good grades. For Perseus, 
the fraction of bad grades exceeds 15 per cent even for the 
observations taken with a reduced frame rate (ObsID 3404), 
and for H1821+643, the observations are heavily affected 
by pileup such that a readout str eak is also se e n. Inst ead, 
we use the nuclear flux quoted in Russell et al.l (|2010al ') for 
H1821-I-643 and in iMerloni fc Heinj (|2007l )" for Perseus, the 
latter of which is an average estimate based on the fluxes 
available in the literature. The nuclear flux of Perseus also 
appears to have varied significantly over the last 10 years, ac- 
cording to the different values published in the literature (see 
Allen et al.l I2OOII: lEvans et al'l I2OO6I : IChurazov et al.l l2003l : 



Donato et al.l 20041 ). To account for the order of magnitude 



variability, we assign an uncertainty of 0.5 in logarithm space 
(log Lx = 43.4 ±0.5). 

Applying a similar method, this time for a 12.7kpc 
X 12.7kpc square region which comfortably contains 



Chandra's PSF, we find that the scatter increases slighty, 
but that our conclusions discussed in Section 8 remain the 
same. We also investigate another method for deriving the 
nuclear fluxes. It consists of fitting a regression to the surface 
brightness profile within the central tens of kpc, for exam- 
ple a power-law, and then attributing the extra number of 
counts in the nuclear region as those associated with the cen- 
tral point source. The observe d additional coun ts can then 
be converted to a nuclear flux (| Allen et al.ll200^ ). However, 
as mentioned in HL2011, we investigated this possibility for 
many of our objects at low redshifts, but this method pre- 
dicted a higher number of counts than that found within the 
nucleus. This is because many of our clusters have surface 
brightness proflles that flatten strongly within the inner re- 
gions, and this extrapolation approach cannot be used to 
estimate the flux of the nucleus. 



4.2 Nuclear luminosities based on spectral model 
fitting 

In the previous section, we obtained the core X-ray fluxes 
based on the observed count rates. This allowed us to 
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Table 4. Spectral modelling for the sources with evidence of non-thermal emission. (1) Cluster name; (2) Model name; (3) Absorbing 
column density at the redshift of the source; (4) Power-law index; (5) Gaussian central energy; (6) Gaussian dispersion; (7) Parameter 
for each model; (8) Unabsorbed nuclear X-ray l u mino sity in the 2 — 10 keV energy band. Galactic absorption was included for all sources, 
and was frozen at the values of iKalberla et a~ I 1I2OO5I '). " Covering fraction of the absorber. ''Temperature and abundance values for the 
MEKAL model. '^Power-law inde x of the first, una b sorbe d power-law (at soft X-rays). ''X-ray nucle ar luminosity for the central galaxy in 
the Perseus cluster, taken from lMerloni fc Heind l|2007l ). "^Values for H1821-I-643 were taken from lRussell et al ] l|2010al V 



Name 


Model 




(1022 cm-^) 


r 


-^gaussian 

(keV) 


^gaussian 

(eV) 


Par am. 




Lx nuc(2 - 10) keV 
(1012 erg s-i ) 


(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


Perseus cluster 


















'^2515'^ 


Cygnus A 


I: phabs 
II: phabs 


[ pow + zphabs(pow+ga)] 
[ mekal + zphabs(pow+ga)] 


^ ' -"-3.0 


1 5+0.5 

1 O+0.5 
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50 
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[3.3 keV; 3.2Zo] 
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RXC J1558. 3-1410 


I: phabs 
II: phabs 


pow + zphabs(pow) ] 
mekal + zphabs(pow) ] 




1.9 
1.9 






[2.2 keV; O.OOZq 




5 0+1-3 


RXC J1524. 3-3154 


I: phabs 
II: phabs 


pow ] 

zphabs(pow) ] 




-^"-'-0-2 
3 4+10 










0.81^2 


4C+55.16 


I: phabs 


pow ] 




1 54+0.07 










12.011;; 
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9 0+0.6 
^■"^-0.5 

2-2-0. 7 






[1.4 keV; O.43Z0 




5 g+O-4 
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MACS J1931. 8-2634 


I: phabs 


zpliabs(pow) ] 


1 s+01 


-1 7+0.2 
' -0.2 










53-21^7 


MACS J0947.2+7623 


I: phabs 
II: phabs 


zpcfabs(pow) ] 
zpcfabs( pow + ga) ] 


r 0+0.8 

0..5_Q g 


1 6+°-^ 

l-°-0.2 
1 4+0.4 


Q+0.1 
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[0.95lS;S'2]" 
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i8oi;; 
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MACS J2046. 0-3430 


I: phabs 
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mekal + zphabs(pow) ] 


ir+.'il 
'-^-10 


1.9 
1.9 






[1.8 keV; 0.39Zo 
[1.8 keV; O.SOZq 
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5-91^;^ 


MACS J0913. 7-1-4056 


I: phabs 
II: phabs 


zpcfabs(pow + ga) ] 
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obtain the fluxes in a systematic way for all clusters, which 
is especially useful for those that have no detectable X-ray 
nucleus. In this case, the bacltground-subtracted count rate 
derived in the previous section could still include some 
thermal contribution. The luminosities obtained should be 
regarded as a conservative estimate of the upper limits of 
the non-thermal contribution. However, 13 of our objects 
show evidence of a central point source associated with the 
BCG, and in some cases, there is evidence of significant 
absorption in the nuclear spectra. The true unabsorbed nu- 
clear X-ray luminosities for these 13 sources might therefore 
be significantly higher than those calculated in the previous 
section, and the evolution we find in Fig. |4] might be even 
steeper. The clusters that have a central X-ray point source 
are Cygnus A, RXC J1524.2-3154, RXC J1558.3-1410, 
4C-h55.16, IVIACS J0547.0-3904, ]V[ACS J1931.8-2634, 
IVIACS J0947-f7623, MACS J2046.0-3430, 
IVIACS J0913.7-h4056, MACS J1411. 34-5212 and 
MACS J1423.8-f2404. Note that Perseus and H1821-)-643 
also have a bright central X-ray point source, but the 



Chandra observations are affected by pileup. As mentioned 
in the previous section , we use the nuclear flux quoted in 
iMerloni fc Heind (l2007l) for Per seus, a nd for H1821-)-643 we 
use the one from iRussell et al] (|2010ah . 

The spectral models adopted for each of the 13 sources, 
as well as the resulting, unabsorbed 2 — 10 l:eV nuclear lu- 
minosities are shown in Table ID Appendix A presents the 
details of the spectral modelling. For each source in Table |4] 
we report the average luminosity across the different mod- 
els fitted. The uncertainty on this averaged luminosity is 
the quadratic sum of the errors where we calculate both 
the upper bound and lower bound uncertainty. However, we 
assume that the errors are symmetric and therefore only 
consider the average between the upper and lower bound 
uncertainty. For the remaining 19 objects that have no de- 
tectable X-ray nucleus, we take the 3cr upper bound of the 
2 — 10 keV luminosities shown in Table |3] (Column 4) as the 
upper limit to the non-thermal emission. Combining these 
19 non detections with the 13 detections of Tabled we il- 
lustrate the new evolution in Fig. \5\ 
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Figure 5. Logarithm of tlie 2 — 10 keV nuclear luminosity as a function of redshift for all BCGs in our sample (in units of erg s"^ ). 
The 13 objects that have nuclear luminosities derived from detailed model fitting are shown with the blue circles, whereas the 19 objects 
that have no detectable X-ray nucleus are shown with the downward pointing arrow (3fr upper limit). The Perseus cluster, Hercules A, 
H1821+643, Cygnus A and 3C 295 arc highlighted with the stars. 



To illustrate the steepness of the evolution, we first per- 
form a simple linear regression in the log-log space such that 
^nucleus oc [1 + z)"" . We stress that the scatter in our plots is 
large and various selection effects may be present (see Sec- 
tion 6) . We therefore cannot derive a precise estimate of how 
these BCGs evolve with time. We only use the regression to 
demonstrate that the AGN in BCGs with X-ray cavities are 
on average brig hter at higher redshift. We use the method 
of iKellvl (|2007l n based on a Bayesian approach that uses 
a Monte Carlo technique to simulate the linear regression 
parameters from their probability distribution given the ob- 
served data. Note that upper limits can be incorporated into 
the fit, but these influence strongly the flt and drag the nor- 
malization downwards. Applying this method, we flnd that 
the slope of the regression is positive to more than a 95 per 
cent confldence level (i.e. ^2cr). 



5 RADIO LUMINOSITIES 

In this section, we focus on the radio properties, which can 
also be used as a tracer of black hole activity. Most of the 
radio emission in AGN arises from relativistic particles spi- 
ralling around magnetic flelds and therefore emitting syn- 
chrotron emission. This type of emission is characterised by 
a steep radio spectrum with a > 1 [Sv oc v^"). The radio 
lobes, observed to correlate with X-ray cavities in clusters 



are dominated by this emission. At higher radio frequencies 
(«GHz), emission from the core starts to dominate and is 
characterized by a flat or inverted spectrum with a ~ 0. 
Here, the core refers to the innermost brightest region of 
the jetted emission. We therefore consider more than one 
frequency in our analysis. 

Nine of our c luster s have already been analysed in detail 
bv lBIrzan et al.l (|2008l ) at radio wavelengths. These authors 
derive flux density values at 327 MHz, 1.4 GHz, 4.5 GHz and 
8.5 GHz. They also derive total radio luminosities by fitting 
a model to the spectral distribution of the source between 10 
MHz and 10 GHz. Th ese values are shown in Table [5] Note 
that, for some clusters. iBirzan et al.l (|2008l ) derive both total 
and partial flux densities, the latter of which only includes 
the contribution of the radio lobes. We consider the total 
contributions. The radio properties of the MACS clusters 
are reported in HL2012. These values are also reported in 
Table [5] For the remaining clusters, we proceed in the same 
manner as in HL2012, and searched for radio emission asso- 
ciated with the central AGNs by using various radio surveys 
publicly available. 

First, we searched for evidence of emission at 74 
MHz with the V LA Low frequency Sky Survey (VLSS; 
ICohen et all |2007l ) at 150 MHz with the TIFR GMRT 
Sky Survey (TGS^I) and at 326 M Hz with the Wes t erbor k 
Northern Sky SurvejH (WENSS; iRengeUnk et~all Il997l ). 
Next, we searched for emission at 843 MHz using the Syd- 



* The LINMIX_ERR.PR0 script is available at the IDL astronomy ^ http://tgss.ncra.tifr.rcs.in 

library, |http :/ / idlastro . g sfc . nasa. gov/ ^ http:/ /www.astron.nl/wow/testcodo.php?survey=l 
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Table 5. Radio pr operties - (1) Name; (2) 74 MHz radio flux density from VLSS; (3) 326 MHz radio flux density from WENSS or 
iBirzan et al] bOOSi) ; (4) 1400 MHz r adio flux density; (5) 5 GHz radio flux density; (6) 28.5 GHz radio flux density from Berkeley- 
lUinois-Maryland Association fBIMA: ICoble et al.ll2007l ): (7) Other radio flux density; (8) Integrated radio lu minosity from 10 MHz t o 
10 GHz; (9) References: (i) 1.4 GHz FIRST; (ii) 1.4 GHz NVSS; (iii) Parkc s MIT-NRAO 4.85 GHz survey (Griflitli fc Wrightj |l993l) ; 
(iv) 5 GHz private communication (M T. Hogan); (v) 5 GHz fromlHines fc Wi l ]s (199^); (vi) 5 GH z from [CavaKnolo et al. (20 1]]); (vii) 
5 GHz fromlGregorv fc CondonI l|l99j ): (viii) EMSS 5 GHz l|Gioia et al.lll990l : Istocke et al.^ll991^ : (ix) 843 MHz SUMSS ijBock et al'l 
19991; iMauch et al.ll2003l) ; x) 15 GHz Arcminute Microkelvin Imager (AMI), private communication (K. Grainge); (xi) 1.4 GHz from 
Ehlert et al.l ll201ll ); (xii) 150 MHz TI FR GMRT Sky Sur vey; (xiii) Radio fluxes at 327 MHz, 1.4 GHz, 4.5 GHz and 8.5 GHz, as well as 
the total radio luminosities derived in lBirzan et al.l }2008l) . "For the central galaxies of MACS J0242.5-2132, MACS J1411.3+5212 (3C 
295), PKS 0745-19, 4C-I-55.16, A2029 and Hercules A, we use the extensive follow-up radio observations available in the NASA/IPAC 
Extragalactic Database (NED) to determine Liadio- ''Abell 2029 and Hercules A have extended morphologies in the different surveys 
due to resolved radio lobes. The fluxes shown are the total integrated ones for all emission associated with the central AGN, including 
the contribution of the radio lobes. 
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0.16 


(ii,ix,xii) 


PKS 0745-19° 


39060±200 




2372+119 


480+27 




26770+1350[O150MHz) 


6.07 


(ii,iii,xii) 


Abell 2204 


1640±290 




57.9+2.9 


25.9+1.3 


8.79+0.13 




1.23 


(i.iv) 


Hercules A"'** 


1038200±300 




47120+2356 


11376+592 






313.7 


(ii) 


Abell 0115 


38810±304 




1362+68 
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0.40 


(xiii) 
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(xiii) 
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3.2 
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< 116 
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0.42 


(ii,viii) 
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1.02 


(ii,vi) 
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0.94 


(i.iv) 
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1.23 


(i) 
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7.2 


(ii) 


MACS J0159.8-0849 
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31.4+ 1.6 
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(i.iii) 
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(ii,ix) 


MACS J0913.7+4056 


< 272 


54.0 + 7.1 


8.3 + 0.5 


1.6 + 0.1 


0.69 + 0.12 


0.80 + 0.04[ei5GHz] 


l.,33 


(i,v,x) 


MACS J1411. 3+6212° 


120270 ± 6022 


61647 + 3082 


22171 + 1109 


7401 + 808 






1025.3 


(i.vii) 


MACS J1423.8+2404 


< 232 




5.2 + 0.4 




1.49 + 0.12 




3.58 


(i) 



ney University MqlonRl o Sky Survey (SUMSS; iBock et all 
ll999l : lMauch et alj[200a ). We determined radio fluxes at 1.4 
GHz using the 1.4 GHz VLA Faint Images of the Radio 
Sky at Twenty-Centimeters survey (FIRST; iBecker et all 
1 1994 ). If FIRST data were not available, we derive the 
fluxes u sing the 1.4 GHz NR AO VLA Sky Survey catalogue 
(NVSS: [Condon et alj[l998l ). Finally, we searched for high 
frequency radio emission associated with the central AGNs 
using the 5 GHz Parke s-MIT-NRAO (PMN) radio survey 
(|Griffith fc Wrightll993h and the compilation of results from 
EINSTEIN Observatory Extended Medium-Sen sitivity Sur- 
vey (EMSS: [Gioia et al.|[l990l : IStocke et al.lll99ll ). Finally, at 
28.5 GHz, we used the Berkeley-Illino is-Maryland Associa- 
tion survey (BIMA: [Coble et al.|[2007h . As in HL2012, if no 
point source was seen within 100 kpc of the BCG, we used 
the 2arms value within the beam area as an upper limit to 
the flux. 

For some of the radio surveys, the errors quoted in the 
catalogue did not account for systematic uncertainties. Sys- 
tematic uncertainti es vary with frequ ency, but are on the 
order of 5 per cent (|Carilli et aLlll99l] '). If a survey did not 



include systematic uncertainties, we computed the total un- 
certainty assuming a 5 per cent systematic error and a 2Gr-mB 
noise level. In this case, errors are derived as the quadratic 
sum of the rms noise level in the map and the systematic 
uncertainty associated with the value. 

Finally, we computed a rough estimate of the total radio 
luminosity between = 10 MHz and V2 = 10000 MHz 
following Eq. \^ where Dl is the luminosity distance to the 
source and 5*^ is the flux density. 

Lradio = 47rl>£ f \s.)dv (2) 

The flux densities at 10 MHz and 10 GHz are extrapolated 
from the other known values based on the assumption that 
the flux density scales as Sv oc v~" . Here, the local spec- 
tral index determined from the two nearest flux density 
data points are used to compute the extrapolated values. 
Using a simple trapezoid rule, we then integrated over the 
10 MHz to 10 GHz range. For MACS J0242. 5-2132, MACS 
J1411.3-F5212, PKS 0745-19, 4C-F55.16, A2029 and Hercules 
A, we also used the extensive follow-up radio observations 
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Figure 6. Left - Logarithm of the total radio luminosity in erg as a function of redshift over the 10 MHz to 10 GHz range. The 
Perseus cluster, Hercules A, H1821+643, Cygnus A and 3C 295 are highlighted with the stars. Right - Same but focused on the radio 
emission seen at 5 GHz. Also shown in the right panel are the errors bars, which have been multiplied by a factor of five for illustrative 
purposes. 



available in the NASA/IPAC Extragalactic Database (NED) 
to determine Lradio- Our results are shown in Fig. |6] Since 
high frequency radio observations are better proxies of the 
core emission, we also show in Fig. [S] the 5GHz radio lumi- 
nosities of our sources as a function of redshift, but stress 
that the scatter in the radio figures is larger than that ob- 
served at X-ray wavelengths (Fig. [S]). We therefore concen- 
trate only on the X-ray evolution of our sources in the fol- 
lowing sections. 



6 RESULTS AND SELECTION EFFECTS 

Fig. E] shows that the X-ray nuclear emission of BCG- hosted 
AGNs is evolving with time such that they are on average 
brighter in the past. These BCGs have been selected such 
that they have known clear AGN-driven X-ray cavities, and 
all lie in cool-core clusters. They therefore represent a subset 
of the BCG population and the evolution seen may not nec- 
essarily apply to all BCGs. In other words, we only consider 
those with strong radio-mode feedback taking place, and 
trace how the radiative evolution of these systems evolves 
with time. 

We stress that there is no clear selection bias against 
the lower right portion of Fig.O i.e. the portion where high- 
redshift BCGs have low nuclear X-ray luminosities. Indeed, 
when selecting our sample, our main criteria is whether a 
massive cluster has X-ray cavities. The detectability of X-ray 
cavities does not depend on the central AGN being switched 
on or not radiatively, since Chandra can resolve the nuclear 
point source from the cavities. If the nuclear X-ray emis- 
sion of 3C 295 was 100 times fainter, we still would have 
included it in our sample since we still would have seen the 
X-ray cavities. Yet, the lower right portion of Fig. [5] remains 



unpopulated, pointing to some form of evolution of the nu- 
clear properties in BCGs with X-ray cavities. In this section, 
we investigate other possible selections effects that may alter 
our results. 

First, we discuss BCG duty cycles in cool-core clus- 
ters. Studies of nearby clusters show that almost all cool- 
core clust ers have X-ray cavitie s (detect i on ra tes > 90 per 
cent, see iDunn fc" FabianI l2006l : iFabianI |2012| ). suggesting 
that BCG-hosted AGN duty cycles are large (> 90 per cent) 
in cool-core clusters, at least in terms of mechanical out- 
flows. A possible selection effect could therefore be that we 
are missing some clusters that currently have no apparent 
X-ray cavities and therefore have not been included in our 
sample due to our selection criteria, but these represent at 
most 10 per cent of the cool-core population and should not 
affect our results significantly if the duty cycles remain high 
at high redshifts. 

An important selection effect is that we could be missing 
systems with X-ray cavities at high redshift, not due to the 
AGN being switched off in terms of mechanical outfiows, 
but rather due to the increasing difficultly in detecting X- 
ray cavities at high redshift combined with a lack of deep 
data for these objects. In HL2012, we found that at least 
19 out of 37 MACS clusters with cool cores (50 per cent) 
seemed to host cavities. These 19 clusters included 13 with 
clear cavities as well as 6 with less well defined cavities and 
only classified as potential cavities (see IIL2012 for more 
details). More importantly, the majority of the remaining 
18 cool-core clusters with no detectable X-ray cavities have 
less than 30 ks Chandra observations. A possible selection 
effect could therefore be that we are missing a significant 
population of MACS clusters with X-ray cavities for which 
the data quality was not sufficient to detect cavities but 
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deeper data for these objects would be needed. While being 
mindful of this possible effect, for now, we base our results on 
the current population of BCGs with known X-ray cavities. 

Third, it is possible that we are missing Lx-ray > 3 x 
10*"* erg s~^ clusters at high redshift not initially included 
in MACS due to the flux limited nature of the survey. New 
surveys based on the Sunyaev-Zel'dovich effect which has 
no surface brightness dimming, are starting to come online. 
These include those based on the South Pole Te lescope (SPT 
ICarlstrom et al ] |201ll : IVanderlinde et al.ll20ld ) and Planck 
(Bartlett et al. 2009, Planck Collaboration 2011a), and will 



provide a wealth of new data in the intermediate to high 
redshift Universe. 

Note that, as we are increasing in redshift, we are prob- 
ing a larger comoving volume. A possible selection effect 
could be that if clusters with radiatively efficient central 
AGN are rare, as we increase in redshift, the probabil- 
ity of encountering one increases. This could explain why 
we detect such objects at high redshift, but can not ex- 
plain the lack of very faint ACN at high redshifts which 
should remain just as common if there was no evolution. 
The evolution we find is also steeper than the expected in- 
crease in comoving volu me with redshift for our cosmol- 
ogy. based on the iRowan -Robinsori (196^) test (see also 
iRowan-Robinson fc Wand |2010| ). so this effect cannot ex- 
plain entirely our results. 

Another possible selection effect could be that we are 
selecting only the brightest clusters at high redshift, since 
these are the ones where X-ray cavities are more easily de- 
tected. However, the top panel of Fig. [7] shows that even 
if we correct the cluster X-ray luminosity, we still see that 
higher-redshift BCCs are more X-ray active than their low- 
redshift counterparts. Here, we use the cluster luminosities 
shown in column 4 of Table [2] In this linear regres- 

sion in the log-log space is consistent with a positive slope 
at a 92 per cent level. 

We also investigate the possibility that the evolution 
seen is not due to the intrinsic evolution of the nuclear lu- 
minosities, but instead caused by the evolving properties 
of the cluste r. For lower r edshif t clusters at 2 < 0.5, both 
HL2012 and iBauer et ahl (120051 ) show that the fraction of 
cool-core clusters to non cool-core clusters remains fairly 
constant up to present day, with almost half of the clusters 
showing evidence of short central cooling times. No evolution 
in the cool-core properties of clusters are seen, at least up 
to z = 0.5 (however, see Samuele et al. 2011 and McDonald 
2011 for evidence of evolution in terms of the optical proper- 
ties of BCGs). At 2 > 0.5, several authors have shown that 
their appears to be a significant lack of very strong cool- 
core clusters (IVikhlinin et al.ll2007l : ISantos et al]l2008l . [201ol : 
lAii en et al.ll201ll ). indicating that there may be some form 
of rapid evolution in the cool-core prop erties of clusters be - 
tween 2 = 0.5 and 2 = 1. However, see ISemler et al] (|2012D 
for recently identified 2 > 0.5 strong cool core clusters based 
on the SPT survey. Our sample consists of clusters within 
2 = 0.6. At least in terms of the cool core properties of the 
clusters, there should therefore be no significant evolution, 
and this should not explain the trend we see in Figs. |4] and 
[S] However, we have computed the cooling luminosities of 
all of our clusters in Table [2] and show in the middle panel 
of Fig.[7]the nuclear X-ray luminosities corrected for cooling 
luminosity. The linear regression remains consistent with a 
positive slope to a 90 per cent. Note however, that the ob- 
served lack of very strong cool-core clusters at 2 > 0.5 could 
also be due to an identification bias (see Section 8.3). 

According to Table [2j our 2 > 0.3 clusters have cooling 
luminosities (Lcooi) 2 to 3 times higher than our 2 < 0.3 
clusters. Larger cooling luminosities implies the need for 
stronger feedback from the central AGN to counterbalance 
the cooling, and may therefore result in a higher accretion 
rate and brighter AGN. The observed evolution in Fig. [S] 
could therefore simply be due to a selection bias where we 
selected on average weaker cool cores at low redshift. We 
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Figure 8. Same as Fig. [S] but we limit our sample to 26 clusters 
instead of 32. We remove the weakest cool core clusters at 2: < 0.3 
(Abell 0085, Cygnus A, RXC J1524.3-3154, RXC J1558.3-1410, 
Hercules A, Abell 0115 and ZwCl 2701) such that the average 
cooling luminosity for our z < 0.3 clusters is the same as the 
one for our z > 0.3 clusters. In other words, we further limit the 
sample so that we are looking at the same population of cool core 
clusters with X-ray cavities across all redshifts. 



test this theory in Fig. [51 where we hmit the sample to 26 
clusters instead of 32 such that the z > 0.3 population of 
clusters has on average the same cooling luminosity as the 
z < 0.3 population. To do this, we simply remove the weak- 
est cool core clusters until the average cooling luminosities 
are within 10 per cent of each other. This results in the re- 
moval of Abell 0085, Cygnus A, RXC J1524.3-3154, RXC 
J1558. 3-1410, Hercules A, Abell 0115 and ZwCl 2701, aU 
of which are the weakest cool core clusters in our sample, 
i.e. the clusters with the smallest cooling luminosities. If we 
limit the sample to the remaining 26 clusters, and plot the 
2 — 10 keV nuclear luminosity as a function of redshift for 
this limited sample, we see from Fig. [8] that the apparent 
evolution is even more pronounced. Here, the slope of a lin- 
ear regression is positive to more than a 96 per cent level. 

Finally, we analyze in the bottom panel of Fig. [7] 
the nuclear X-ray luminosities corrected for cavity energy, 
where we estimate that the energy stored within each 
cavity is given as the multiplication between the thermal 
pressure of the surrounding hot gas and the cavity vol- 
ume (PVcavitics). The cavity energetics sho wn in Fig. [7| 
were taken from the litera ture (see Table [2l Birzan c t al. 
2004 : Dumi ct al. 2005; Raffertv et"al]|2006l;l&un n & Fabian 
20061 . l2008l : iHlavacek-Larrondo et akl l201ll . and HL2012), 
but RXC J1524.3-3154 or RXC J1558.3-1410 have not 
been included in the plot since the cavity energetics are not 
yet available. For the objects where multiple cycles of out- 
bursts can be seen (e.g. Perseus), we only consider the inner 
X-ray cavities. Fig. [7] shows that when correcting for cavity 
energy, the scatter is larger but the slope remains positive, 
indicating that the evolution seen is mostly due to the evo- 



lution of the radiative properties of the BCGs. We further 
discuss this plot in Section 8.1.3. 

Note that, we have used the cavity enthalpy as a proxy 
of the outflow energetics, and not the cavity power (Pcavity ~ 
4PVcavitios/iagc) as the literature often shows different es- 
timates of a cavity age (e.g. sound crossing time and/or 
buoyancy rise time), and it remains unclear which one is 
the most accurate. To avoid the added scatter of cavity age, 
we have chosen to illustrate Fig.[7]using simply the enthalpy 
(f Vcavitios). Note also, Perseus is the only source that clearly 
shows variability of the nuclear emission over yearly time 
scales, and we include an error bar for this source that takes 
this variability into account. It is possible that other BCGs 
also exhibit such variability, but it is not clear if this is the 
case for aU BCGs. 4C-F55.16 and PKS0745-19 are two ex- 
ample tha t show no significant variation o n 5-10 year time 
scales, see iHlavacek-Larrondo et all l|201lh and Sanders et 
al. in preparation. 

In summary, selection effects are important, and while 
being mindful of them, we base our discussion in the fol- 
lowing sections only on the currently known population of 
BCGs that lie in massive clusters of galaxies with known 
X-ray cavities. 



7 HIGH-REDSHIFT CLUSTERS OF 
GALAXIES 

There are a number of recent studies that have analysed 
individual 2: ~ 1 clusters, and found evidence for both 
AGN activity and short central cooli ng times. These in - 
clude WARP J1415. 1-^3612 at z = 1.03 dSantos et al.ll2012l ). 
XMMUJ0044.0-2033 at z = 1.579 (ISaritos et al.ll201l( ). PKS 
1229-021 at z = 1.04 llRusseU et al.| 20121 ) "and 3C186 at 
z = 1.067 (ISiemiginowska et al.l 201C ). Although none of 
these studies have provided clear evidence for the existence 
of X-ray cavities being carved out by the radio lobes, most 
likely due to the limited resolution, the X-ray luminosities 
of these clusters are on the order of ^lO*** erg s~^ and all of 



these show at least tentative evidence of a cool core. These 
clusters are therefore similar to our sample of objects, and if 
BCG-hosted AGN duty cycles remain high in high-redshift 
cool-core clusters, most should harbour X-ray cavities. We 
therefore plot these in Fig. [9] Note that we did not include 
these in our original sample since they do not show evidence 
of X- ray cavities. F i g.[9]d oes not show XMMUJ0044.0-2033 
since ISantos et al.l (|201ll ) do not derive flux estimates for 
any non-thermal contribution of a central AGN, i.e. they 
only focus on deriving fluxes for the thermal cluster emis- 
sion. There are also no Chandra observations available for 
this source in order to isolate the nuclear emission of the 
BCG. 

As for WARP J1415.1-H3612. ISantos et~all (|2012l ) anal- 
ysed very deep Chandra observations of this source (278 ks), 
but did not analyse the X-ray emission from the central 
AGN. We have reprocessed the Chandra observations (Obs 
ID 12255, 12256, 13118 and 13119), and searched for evi- 
dence of non-thermal emission associated with the central 
AGN using the method outlined in Section 4.2. In this case, 
we considered a region centred on the BCG within a radius of 
1" as our nuclear region and a surrounding annulus located 
within r — 2" — 3" as our background. Applying C-statistics 
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Figure 9. Nuclear X-ray luminosities as a function of redshift. The black filled circles illustrate the data points from Fig. [5] The downward 
pointing arrows represent upper limits to the nuclear luminosities. The high-redshift sources discussed in Section 7 are shown with the 
blue points. The black dashed curve shows observed evolution seen in star formation efficiency (L* oc (1 + z)"^). The right Y-axis shows 
the Eddington ratio of the nuclear X-ray luminosity assuming a 10^" Mq black hole and the striped region shows the location where a 
black hole reaches quiescence (inucleus'SlO^^LEtjtj)- 



to account for the low number of counts, we find that both an 
absorbed (Galactic) power-law and mekal model provide an 
equally good fit to tlie data. Tliere is no clear evidence from 
the spectrum that the source contains non-thermal emission. 
Fitting a combined mekal -f power-law model to the data 
does not improve the fit, and even when freezing most pa- 
rameters, the fit remains uncontrained. We therefore derive 
an upper limit to the non-thermal emission using the same 
method outlined in Section 4.1 where we convert a count 
rate into a flux. Throughout the analysis, we have assumed 
the same gal actic absorption and central location for our 
regions as in ISantos et all (|2012l . ; iVh = 1-05 x 10^°cm^, 
RA=14:15:11.08 and DEC=-h36:12:03.1). We derive a So- 
upper limit of 1.8 X 10*^ erg s~^ for the nuclear emission of 
WARP J1415.1-f3612. This result is also shown in Fig. [5] 

Fig. [§] also includes the powerful quasar and 
hyper- luminous infr ar ed galaxy I RAS F15307-f 3252 
(|lwasawa et al.l l2005l : iFabian et al.l 1 19961 ). which is a 
Compton-thick AGN embedded in a luminous cluster 
(ix-ray, cluster ~ 10** erg s~^ ). The estimated 2 — 10 keV 
luminosity is on the order of 10*^ erg s~^ . Although there 
is yet no known evidence indicating the presence of a 
cool core or X-ray cavities, we include it in Fig. [9] for 
comparison. Finally, we include the z — 0.596 extreme 
cluster SPT-CLJ2344-4243 in Fig. M also known as the 
Phoenix cluster. This cluster is highly X-ray luminous 
(Lx-ray, cluster ~ ]^q45-46 ^^.g ^-1 -j ^ harbours a stroug cool 



core as well as a massive starburst at its core (74OM0yr ^; 
iMcDonald et all I2OI2I) . No X-ray cavities have yet been 
reported in this object and we therefore only include it in 
Fig. [9] for illustrative purposes. 

Interestingly, Fig.[9]shows that if we were to extrapolate 
the apparent evolution seen in our sample of BCG-hosted 
AGN out to 2 = 1, some of the currently known high-redshift 
clusters with evidence of AGN activity and short central 
cooling times lie within this extrapolation. At least in terms 
of the radiatively efficient high-redshift sources, these ob- 
jects may be the less evolved counterparts of nearby massive 
cool-core clusters with powerful AGN-driven outfiows. Note 
that WARP J1415. 1-1-3612 may simply be absorbed and 
therefore undetectable at X-ray wavelengths with Chandra. 
It therefore does not necessarily represent a distinct popu- 
lation of objects. 



8 DISCUSSION 

8.1 Implications for the transition between 
quasar-mode and radio-mode feedback 

In HL2012, we analysed the properties of X-ray cavities in 
the M ACS sample and com pared these to low-redshift sam- 
ple of iRaffertv et all ()2006l ). Interestingly, we found no sig- 
nificant evolution of the X-ray cavity properties with red- 
shift. We illustrate this in terms of the cavity energetics in 
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Figure 10. Left: Cavity enthalpy as a function of redshift for the sources in our sample (see Table[2]|. Middle: Fraction of BCGs in our 
sample with detected X-ray nuclei as a function of r e dshift (blac k). Also shown are the fraction of galaxy clusters with X-ray cavities 
based on the low-redshift sample of iDunn fc FabianI ||2006| . l20od ; 30 per cent) and the z > 0.3 sample of HL2012 (rj 20 - 30 per 
cent). For the latter, we show the fraction in terms of the clear and potential cavities (light blue) and the fraction in terms of only the 
clear cavities (purple). Right: Histogram of the number of sources in our sample with a detected X-ray nucleus. Note that, our sample 
only includes BCGs that lie in massive clusters and that have known X-ray cavities. 



the left panel of Fig. IIOI where we plot the cavity enthalpy 
as a function of redshift for the sample studied here. We 
stress that our sample only includes massive clusters with 
a 0.1 — 2.4 keV luminosity greater than 3 x 10*^ erg s~^ . 
Fig. [To] shows that there is no significant evolution in the 
cavity energetics (less than an order of magnitude), at least 
in terms of the cavities found in massive cool core clusters. 

On the other hand, Fig.[9]shows that BCGs with X-ray 
cavities appear to be more radiatively efficient at high red- 
shifts in terms of their X-ray luminosities, suggesting that we 
are seeing some form of evolution of how AGN feedback op- 
erates in BCGs with X-ray cavities ov er the z = 0—0.6 range. 
Such an evolution was first noted bv lHall et all (| 19971 ') . who 
analyzed a small subset of 3 quasar-like BCGs and found ev- 
idence for evolution. Here, we consider a significantly larger 
sample of BCGs and analyze the implications of our results. 
We investigate two possible scenarios that explain our re- 
sults, and dicuss them in the following two sections. 



8.1.1 Evolution in terms of the general increase in 
luminosity with redshift 

We first consider the possibility that we are seeing a general 
increase in the X-ray luminosity of BCGs with increasing 
redshift (see top panel of Fig. Ilip . in the sense that BCG- 
hosted black holes are rapidly evolving such that their 2 — 
10 keV X-ray luminosity has increased with redshift by an 
order of magnitude over the last 5 Gyrs. 

Our sample of BCGs only includes those in the most 
massive clusters of galaxies and that have clear AGN-driven 
outflows in the form of X-ray cavities. We are therefore only 
considering a subset of the BCG population. We also stress 
that we are only probing the 0.5 — 7 keV energy range as 
seen from Chandra. It is not clear if the evolution seen can 
be applied to all wavelengths. 

Nevertheless, the suggestive evolution seen in Fig. [5] 
is consistent with being steep. In comparison, the star 
formation efficiency in galaxies evolves approximately as 
Li, (X (1 -f z)'^ over the redshift range 2: ~ — 1 (e.g. 



Ruiopakarn et al.ll2010l : iLillv et al.|[l99g : iMadau et ai]ll996l : 
Hopkinj|2004l ). A similar evolution is seen in terms of the 
number density of quasars which peaks at a reds hift of z ~ 2 



and then rapidly declines to the present day facd a et al 



J.003 ; Hasingcr ct al. 2 0051; iHopkins et"ai] [2007 ; Aird et al 
20101 ). The exact shape of density function depends on the 



luminosity bin considered (higher luminosity AGN peak at 
higher redshift), as well as the wavelength considered, but 
the decline from z « 1 to the present day follows roughly a 
(H-2:)^ scaling at all wavelengths (e.g. [Hopkins et al.|[2007l ). 
Fig. [U] highlights our data points as well as the typically ob- 
served (1 + z)'^ scaling with the dot-dashed line. If we cal- 
culate the average nuclear luminosity for our 0.3 < z < 0.6 
BCGs, we find that it is at least 10 times higher than that of 
the BCGs in the < z < 0.3 redshift range, yet the star for- 
mation efficiency increases only by a factor of 2 — 3 over these 
redshift bins. BCGs would therefore appear to be changing 
more rapidly than the general population of AGNs and the 
decline of star formation activity in galaxies. It is possible 
that the unique environments of BCGs, and especially those 
in our sample which lie in the most extreme cool-core clus- 
ters, cause their central AGN to shut down radiatively more 
rapidly than others, perhaps due to higher merger rates at 
the centres of clusters. 

Quasar-mode feedback and its evolution have also been 



implemented in various simulation works (iDi Matteo et al 



2005 



2005 



Croton et al. 2006; Bower et all 200i 



Springel et al 



Siiacki fc Springeli .20061 ; [Merloni fc Hein j |2008|).^ 



partic ular, both Croton et al.[ (|2006l ) and Merloni fc Heinj 
(|2008h predict that the black-hole accretion rate density (re- 
lated to luminosity density assuming a certain efficiency) for 
black holes operating in quasar-mode feedback declines at 
most by a factor of ten since z ~ 2. Furthermore, their sim- 
ulations predict that the accretion rate density for black 
holes operating in radio-mode feedback (or kinetic-mode 
feedback) remains roughly the same since 2: « 2. Although 
both these pictures agree well with our BCGs in terms of 
radio-mode feedback since we find no evidence for evolution 
in the outflow mechanical properties in HL2012, the evolu- 
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Figure 11. Two possible scenarios to explain the evolution of the nuclear X-ray luminosities in our sample of BCGs with X-ray cavities 
as a function of redshift (redshift increases from left to right). The orange circles illustrate AGN that have a detectable X-ray nucleus, 
whereas the black colors illustrate the BCGs with no detectable X-ray nucleus. The top panel shows the first evolution scenario, such 
that BCGs with X-ray cavities become brighter with increasing redshift (from left to right). As they become brighter, they become more 
easily detected, hence the increase in the number of orange points. The bottom panel shows the second scenario, in which the fraction 
of radiatively efficient sources increases with redshift, from 30 per cent at 2 ~ 0.1 (left) to some 60-70 per cent at z 0.6 (right). 



tion of quasar-mode feedback may be significantly steeper 
for BCGs with X-ray cavities than predicted by these simu- 
lations, i.e. a factor of 10 since z ~ 0.5 as opposed to a factor 
of 10 since 2 ~ 2. At least in terms of BCGs in massive clus- 
ters with X-ray cavities, our results suggest that they evolve 
at rates much faster from quasar-mode to radio-mode feed- 
back than predicted by simulations. 

8.1.2 Evolution in terms of the number of sources 
switched-on radiatively 

We also investigate a second possibility as to why BCGs 
with X-ray cavities appear to be more radiatively efficient 
at high redshifts. Our sample only contains a small number 
of sources, especially at high redshifts, and could therefore 
be affected by small number statistics as well as selection 
effects as seen in Section 6. Yet, we find that the high redshift 
sources are brighter, and that the majority have a detectable 
X-ray nucleus. 

If there was no evolution and assuming that we have 
not been biased towards only radiatively efficient sources at 
high redshift, we should have detected the same fraction of 
sources with bright, radiatively efficient nuclei and therefore 
easily detectable at all redshifts. Yet, we find that the ma- 
jority of high redshift sources have a bright and detectable 
central AGN, whereas the majority of the low redshift ones, 
even with very deep observations, have no detectable central 
X-ray nucleus. This points to evolution, and we investigate 
the possibility that the evolution seen is simply be due to 
the increasing fraction of sources switched-on radiatively, as 
opposed to a general increase in brightness for the BCG 
population with X-ray cavities with increasing redshift. In 
other words, there would be two states, one radiatively effi- 
cient and one inefficient and as time progresses, BCGs would 
jump from the efficient to the inefficient state, making them 
largely undetectable at X-ray wavelengths. 

We illustrate this second possibility in the middle and 



right panels of Fig. [10] where we plot the fraction of BCGs in 
our sample that have a detectable X-ray nucleus as a func- 
tion of redshift. Here, we have chosen the redshift bins such 
that each contains 6 or 7 BCGs in them. Since our sample 
only contains 32 sources, the plot is affected by small num- 
ber statistics and varies depending on our definition of the 
redshift bins. However, in all cases, we see that the fraction 
roughly doubles over the last 5 Gyrs, and on average varies 
from 30 per cent at 2; ~ 0.1 to 60 per cent at 2 ~ 0.6. This 
result remains even when we consider the limited sample 
of Fig. [S] which only considers the 26 clusters with similar 
cooling luminosities instead of the original 32 clusters in our 
sample. 

The middle panel of Fig. [TO] also shows the fraction of 
clusters with X-ray cavities in two redshift bins (see the 
coloured points) , highlighting the lack of strong evolution in 
the fraction of cluste rs with cavities. Here, w e use the sam- 
ple of 71 clusters in iDunn fc FabianI (|2008h . which consist 
of a complete sample of cool-core and non cool-core clus- 
ters, to illustrate the average fraction of clusters with cav- 
ities at low redshift (red point), and compare this to the 
MACS clusters at high redshift (blue points). This figure 
therefore shows that the number of clusters with X-ray cav- 
ities remains roughly the same, yet those with cavities and 
detectable X-ray nuclei varies strongly with redshift. A pos- 
sible scenario could therefore be that we are not seeing a 
gradual increase in X-ray luminosity of all of our sources, 
but rather an increase in the fraction of sources that are in 
a radiatively efficient state (see lower panel of Fig. Illjl . 

8.1.3 Brief note on radiatively efficient BCGs with X-ray 
cavities 

The bottom panel of Fig. [7] illustrates the ratio between 
the nuclear X-ray luminosity and mechanical energy of the 
AGN-driven outfiows in our sample as a function of redshift. 
An interesting result is that if we only consider the sources 
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with detected AGN in our sample, the ratio between the 
radiative and mechanical output seems to remain fairly con- 
stant with redshift, apart for the only quasar in our sample, 
H1821-I-643. This might be indicating that once the nucleus 
is switched on radiatively, for instance once it reaches a ra- 
diatively efficient state, the power that emerges in the form 
of mechanical jets scales roughly proportionally with the ra- 
diative power (i.e. how much mass the black hole is accret- 
ing), regardless of redshift. This is discussed in more detail 
in Russell et al. submitted. This result might be indicating 
a fundamental property of accretion physics. 



8.2 Implications for the accretion modes in 
BCG-hosted black holes: analogy with 
black-hole binaries 

We now discuss the possibility that we are seeing a state 
transition of BCG-hosted black holes, operating from a more 
radiatively efficient state at high reds hift to a radiatively in - 
efScient state at low redshift (see also lChurazov et al.ll2005l ). 
This interpretation can be applied to both scenarios pro- 
posed in Fig. 111! 

The right hand axis of Fig. [9] shows the nuclear 2 — 
10 keV luminosity of each object in terms of the Eddington 
ratio. To compute the Eddington ratio, we have assumed a 
typical black-hole mass of 10"' M0 for our sample of BCGs 
since they lie in the most luminous clusters of galaxies, i.e. 
BCG mass and therefore SMBH mass scales with cluster 
X-ray luminosity. It is therefore reasonable to assume that 
these BCGs will host some of the most massive black holes 



with regards to other BCGs and wil l therefore have typical 
average masses of 10^°Ma (se e also iMcConnell et al.ll201ll : 
iHlavacek-Larrondo et aLll2012l ). 

According to Fig.O low-redshift BCGs with X-ray cavi- 
ties are extremely radiatively inefficient and have typical Ed- 
dington ratios of 10"^ — 10~^LEdd. On the other hand, our 
results predict that the equivalent population of BCGs at 
z = 1 will have typical Eddington ratios of 10~^ — 10~''LEdd- 
In analogy with black-hole binaries, such a strong Edding- 
ton ratio change corresponds to a state transition between 
the canonical "low/hard" state and extremely radiatively 
inefficient "quiescent" state. 

Black-hole binaries are important tools in X-ray as- 
tronomy, often referred to as X-ray binaries and when ra- 
dio loud, microquasars. In particular, observations show 
that these objects transit between various accretion states 
on time scales of days/years, and such short time 
scales allow detailed studies of how black holes evolve 
wi th time through monitoring campaigns (see reviews 
bv iMcClintock fc Remillardll2006l : iRemillard fc McClintockl 
I2OO6I ). Black-hole binarv svstems appear to typically tran- 
sit between four major states (excluding the intermediate 
states): a very high state also known as the "steep power- 
law" state which is often accompanied by quasi-periodic 
oscillations (QPOs), a thermal state also known as the 
"high/soft" state characterized by an efficiently accreting 
black hole (^110~^LEdd) in which the jet is quenched, a hard 
state known as the "low/hard" state where the black hole 
is accreting at low rates (^lO^'^LEdd) and driving a pow- 
erful jet, and finally the quiescent state where a black hole 
is acrreting at rates orders of magnitude bellow the above 
rates. The luminosities associated with the quiescent state 



of black-hole binar ies typically range between L x-rav = 
^q30.5-33.5 gj.g g-i (|Remillard fc McClintoc3l2006l ). equiva- 
lent to Eddington ratios in the range of 10~* — 10~^LEdd for 
stellar mass black holes. Some studies have suggested that 
the quiescent state is simply an extensi on of the lo w /hard 
state with even lower accretion rates. iGafio et al.l (|2006D 
showed that the fundamental plane of black-hole activity, an 
established relation between the mass. X-ray and radio lu- 
minosity of a black hole in the low/hard state, extends down 
to quiescence for black-hole binaries. Advection-dominated 
accretion flows (ADAFs) which have been able to explain 
the behaviour of black holes accreting at rates below a few 
per cent of Eddington have also bee n able to explain the be- 
havior of quiescent black holes ( e.g. iNaravan fc McClintockl 
I2OO8I : iNaravan et al.lll996l . ll997l '). including Sgr A*, the qui- 
escent black hole a t the center of our galaxy that is accreting 
at M= 10"^MEdd (|Yuan et al.ll2003l ). 

Based on our knowledge of state transitions in black- 
hole binaries and the strong evolution seen in terms of Ed- 
dington ratios in BCGs (Fig. we propose that BCG- 
hosted black holes with X-ray cavities have migrated from 
the canonical "low/hard" state to the "quiescent" state over 
the last 5 Gyrs. This would explain why the AGN in BCGs 
appear to be so radiatively inefficient at the present day. 
Quiescent states in black-hole binaries are lon g-lived; the 
objects spend most of their time in this state (|Done et al.l 
120071 : INaravan fc McCfintockl I2OO8I . 90 per cent). If BCG- 
hosted black holes with X-ray cavities are transiting into 
this state, then the black holes in these BCGs could pro- 
vide a long-lived solution to the cooling-flow pr oblem if jets 
are m aintained throughout this state (see also iDunn et al.l 
|2010D . 

The transition between the thermal state to the 
low/hard state in black-hole binaries, as opposed to the 
transition between the low/hard state and the high/soft 
state, typically occurs at lower ac cretion rates, 0.02 — 
O.OSLEdd as opposed to 0.1 — 0.2LEdH |Mivamoto et al.|[l995l: 
Maccarone fc Coppi|[200^ : IZdziarskT et al. I I2OO4I : iDone et al.l 



200'if ). At least in terms of BCGs in strong radio-mode feed- 



back with clear X-ray cavities, our results suggest that a 
significant population may host powerful quasars at their 
centres at 2 ~ 1, and some of the currently known z ~ 1 
cool-core clusters fit in this category. These are accreting 
at rates near 10~^LEdd, which places them at the limit of 
the transition state between the high/soft and low/hard 
state analogous to black-hole binaries. This would have oc- 
curred in the last 7 to 8 Gyrs, while black-hole binaries take 
ro ughly 1 year to transit between these states (e.g. G1655- 
10 iMcClintock fc Remillardll2006l : IRemillard fc McClintockl 
l2006h . Scaling the black-hole binary time scale from a 
1 - 5Mq to a 10"'Mq black hole fits well with the observed 
time scale for BCGs. 

We further note that some black-hole binaries in qui- 
escent states exhibit chaotic variations in their luminosities 
(i.e. "fiaring"), while others can r emain stable for decades 
(e.g. A06200-00. ICallo et al.ll2007l ). V404 Cyg is one exam- 
ple and consists of a black-hole binary system that has a 
typical X-ray luminosity of 10^^ erg s~^ in quiescence but 
its lum inosity is known to v ary by a factor of a few within 
hours (|Wagner et all 1 1994 iKong et all |2002| : iHvnes et al.] 
120031 ). Another example is Sgr A*, where the flare fluxes 
can be a factor of a few to ten times those of the quies- 



© 2012 RAS, MNRAS 000, ??-[T8] 



18 J. Hlavacek-Larrondo et al. 



cent state and appear on minut e time scales (|Baganoff et al.l 
l200ll : |Porau"et" et al ] |2003l .[ 20081 ) . The chaotic behaviour seen 
in some quiescent black-hole binaries could be analogous to 
what is seen in Perseus, where the luminosity of the cen- 
tral AGN has varied by a fac tor of a few in the past decade 
(see different flux values by lAllen et al.l I2OOII : lEvans et alJ 
I2OO6I : IChurazov et al.ll200a ). Careful X-ray monitoring of 
this source would be needed to disentangle the variations 
in the nucleus emission. 



8.3 Implications for cluster surveys 

Finally, we mention the implications our results have con- 
cerning cluster surveys and the identification of high-redshift 
cool-core clusters. Based on the results for PKS 1229-021, 
a z ~ 1 quasar e mbed ded in a cool-core cluster of galax- 
ies, iRusseU et all \2Qli ) suggested that high-redshift cluster 
surveys could be missing a significant fraction of cool-core 
clus ters if the central BCG has a quasar-like luminosity (see 
also lAllen et all l201lh . 

Our results suggest that the currently known popula- 
tion of BCGs with X-ray cavities may be evolving with 
cosmic time, such that many may have io*^~*^ erg s^^ 
2 — 10 keV luminosities at 2 = 1. These are quasar-like lumi- 
nosities, implying that a significant fraction of BCG-hosted 
AGN at z = 1 could harbour a quasar in their centres and 
therefore outshine the majority of the cluster X-ray emis- 
sion. Our results therefore suggest that when cluster surveys 
are conducted at high redshift based on moderate resolu- 
tion X-ray imaging (e.g. ROSAT Position Sensitive Propor- 
tional Counters), a significant population of cool-core clus- 
ters might be missed since they would appear as point-like 
AGN and not as extended X-ray objects. This identifica- 
tion bias could explain, at least in part, the apparent lack 
of very strong cool-core clusters at high redshi ft (z > 0.5; 



Vikhlinin et al.l 



20T3). 3C 186 



20071: ISantos et al.|[2008,. 20ld: lAUen et al I 



Siemiginowska et al.|[2010l ') and PKS 1229- 



021 are two such examples of bright quasars embedded in 
cool-core clusters at z ~ 1 where the X-ray emission from 
the BCG dominates over the cluster emission. 

Furthermore, at optical wavelengths, BCGs often ex- 
hibit Seyfe rt-like line emission w ithout being X-ray point 
sources (see lCrawford et al.lll999l . for a detailed look at opti- 
cal properties of BCGs). The X-ray emission associated with 
the Seyfert-type BCG could therefore appear point-like in 
cluster surveys, yet originate from the hot thermal gas of 
the intracluster medium. Existing optical classifications of 
line-emitting BCGs can therefore cause cluster surveys to 
miss cool-core clusters, especially at high redshift. BCGs in 
cool core clusters may also be becoming increasingly blue 
at high redshift, complicating even more the search for such 
objects if the searches are based on the red sequence. 

In summary, several important effects are present here. 
The existence of a significant population of high redshift 
clusters with quasars at their centres therefore complicates 
the search for such objects at high redshifts, at least for 
many of the techniques currently used to find these objects. 
The current population of high redshift cool core clusters 
may therefore be undersampled and t he lack of strong cool 
core clusters observed from surveys (I Vikhlinin et al.l l2007l: 



ISantos et al] |2008| . |201G| : lAUen et all boilh mav simply be 
due, at least in part, to selection effects. 



9 CONCLUDING REMARKS 

We have investigated the evolution of the nuclear 2 — 10 keV 
luminosity in 32 BCGs that lie in highly X-ray lumi- 
nous clusters of galaxies (Lx-ray (0.1 — 2.4 keV) > 3 x 
10''* erg s~^ ) and that have known X-ray cavities. We stress 
that we are only considering BCGs in which strong radio- 
mode feedback is taking place in the form of AGN-driven X- 
ray cavities, and therefore only considering a sub-population 
of BCGs. We further stress that the detectability of X-ray 
cavities most likely depends on the redshift of the sources. To 
limit the possible selection effects, we have applied stringent 
selection criteria to our sample, and only consider highly 
X-ray luminous clusters where X-ray cavities can be more 
easily identified. Applying these criteria, we find evidence 
for evolution such that the average nuclear luminosity in 
these BCGs has increased by at least a factor of 10 from 
< 2: < 0.3 to 0.3 < z < 0.6. If we further limit the sample 
to clusters that have similar cooling luminosities accross all 
redshifts (26 in total), we find that the scatter decreases and 
the increase in brightest with redshift remains the same. 

Mindful of potential biases in our sample, we propose 
that the central AGN of currently known X-ray luminous 
clusters with X-ray cavities are steadily becoming fainter 
over time, or more likely, that the fraction of BCGs with ra- 
diatively efficient nuclei is decreasing with time from roughly 
60 per cent to 30 per cent over the last 5 Gyrs. In analogy 
with black-hole binaries and based on the observed change 
in the Eddington ratios of our sources, we further propose 
that the evolving AGN population in BCGs may be tran- 
siting from a canonical low/hard state analogous to that of 
X-ray binaries to a quiescent state over the last 5 Gyrs. Our 
results also suggest that a significant fraction of BCGs in 
2 ~ 1 clusters may host quasars at their centres, potentially 
complicating the search for such clusters at high redshift. 

New surveys based on the Sunyaev-Zel'dovich effect 
which has no surface brightness dimming (SPT, Planck) 
will soon be online and will provide a wealth of new data 
in the intermediate to high redshift Universe. These, cou- 
pled with follow-up Chandra observations may significantly 
enhance the current sample of clusters with known X-ray 
cavities, thus providing more insight into the potential evo- 
lution observed in this study. Extending the analysis to the 
general population of BCGs, and examining the evolution 
of the X-ray nuclear luminosity across the entire population 
may also provide a better understanding of the evolution 
observed here. 
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APPENDIX A: DETAILED X-RAY SPECTRAL 
MODELLING 

In this section, we present the different nuclear models 
adopted for each of the 11 objects in our sample that show 
evidence of an X-ray point source in their Chandra images. 
These exclude Perseus and H1821-I-643 since we use the nu- 
clear fluxed quoted the literature for these sources. For each 
of the 11 objects, we extract a 0.5 — 7 keV spectrum within 
a 1" circular region centred on the X-ray point source. We 
then take a surrounding annulus located within 2" and 3" as 
the background. All of our fits include Ga l actic absorption 
which we keep frozen at the lKalberla et al] (|2005l ) value. We 
use C-statistics to account for the low number of counts. 



RXC J1524.3-3154 - The background-subtracted nu- 
clear spectrum of this source is very noisy, and requires 
several parameters to be frozen in order to constrain the 
fit. If we fit a simple absorbed (Galactic) mekal model to 
the spectrum, and freeze the abundance at the value ob- 
tained from a mekal model applied to the 2" — 3" sur- 
rounding annulus, we obtain an abnormally large temper- 
ature (I3I7* keV). We therefore investigate other models 
that can explain the observed spectrum, and consider both 
a power-law model without (Model I) and with (Model II) 
internal absorption. The results are shown in Table |4l 

RXC J1558. 3-1410 - The background-subtracted nu- 
clear spectrum of this source cannot be fitted with an ab- 
sorbed (Galactic) mekal model. The fit is largely uncon- 
strained, even if the abundance is kept frozen. If we fit a 
simple absorbed (Galactic) power-law, the power-law index 
tends towards a negative value. This is due to two com- 
ponents being present in the spectrum, one at soft X-rays 
(~ 1 keV) and another at hard X-rays (« 4 keV). Even an 
absorbed (Galactic -|- internal) power-law model cannot fit 
the soft X-rays properly, and the internal absorption tends 
towards a null value. Complex absorbers do not provide good 
fits either. However, fitting a double power-law, where one of 
the power-laws is internally absorbed (zphabs), provides a 
good fit if we keep the photon index of the second, internally 
absorbed power-law frozen at a value of 1.9 (Model I). We 
also investigate the possibility that the emission seen at soft 
X-rays is of thermal origin. If we fit an absorbed (Galac- 
tic) mekal + POWER-LAW model to the data and let the 
temperature, abundance, normalization parameters, as well 
as photon index value vary, the fit is largely unconstrained. 
To provide a better constraint on the model, we adopt the 
method used in HL2011 to derive the non-thermal fiux. In 
this case, the surrounding 2" — 3" annulus is not used as the 
background. Instead, it is used to estimate the properties of 
the cluster thermal component which are then extrapolated 
to the inner 1" circular region. The background is chosen as a 
region located far from cluster emission. For the surrounding 
2" — 3" annulus, we fit an absorbed (Galactic) mekal model 
to the data and find the best-fitting temperature, abundance 
and normalization parameters. The extracted temperature 
is then extrapolated down to the 1" circular region assum - 
ing that T = ar*" where h « 0.3 (| Voigt fc FabianI 120041 '). 
The abundance is not expected to vary significantly from 
r = 1" to r = 2" — 3". Using this abundance and extrapo- 
lated temperature, we first fit an absorbed (Galactic) mekal 
+ POW model to the nuclear spectrum. We let both normal- 
ization parameters and photon index free to vary. We also 
add the constraint that the normalization of the thermal 
component is not allowed to be less than the normalization 
obtained in the fit for the annulus (r = 2" — 3"), scaled 
for the same pixel number, since the density is expected 
to increase with decreasing radius. In this case, the power- 
law index tends towards a negative value. We therefore add 
an internal absorption at the redshift of the source, while 
keeping the power-law index frozen at a value of 1.9 to help 
constrain the fit (Model II). 

4C-I-55.16 - Fitting a simple absorbed (Galactic) mekal 
model to the data provides a good fit, but the resulting tem- 
perature is abnormally high (« 11 keV). This might indicate 
shock heating in the vicinity of the central regions, but we 
also investigate the possibility that the emission is of non- 
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thermal origin since the temperature jump is significant. If 
we fit an absorbed power-law model to the data, where the 
absorption accounts for both Galactic and internal absorp- 
tion at the redshift of the source, we find that the internal 
absorption tends towards a null value. A simple power-law, 
where the absorption accounts only for Galactic absorption 
provides a good fit and we show the results in Table|4] (Model 
I). Fitting more complex absorbers (zpcfabs, pwab) does 
not provide a good fit to the data. 

MACS J0547.0-3904 - The background-subtracted nu- 
clear spectrum is dominated by emission seen around 1 keV. 
We begin by fitting a simple absorbed power-law to the 
background-subtracted spectrum, where the absorption ac- 
counts only for Galactic absorption (Model I) . The resulting 
parameters and unabsorbed nuclear 2—10 keV luminosity 
are shown in Table |4] Although the model provides a rea- 
sonable fit, it underestimates the hard X-rays (> 3 keV). 
Adding an additional absorption at the redshift of the source 
provides a better fit, but the resulting power-law index is 
« 5, which is unlikely. Instead, we interpret the emission 
seen around 1 keV as thermal emission. Fitting a simple 
MEKAL model to the data does not provide a good fit, and 
significantly underestimates the hard X-rays (> 2 keV). 
There is clearly a non-thermal contribution that dominates 
beyond 2 keV. If we fit an absorbed (Galactic) mekal -|- 
POWER-LAW model to the data, and let the temperature, 
abundance, normalization parameters, as well as photon in- 
dex value vary, the fit is largely unconstrained. We therefore 
adopt the same technique as in RXC J1558.3— 1410 to fit the 
data. In this case, the surrounding 2" — 3" annulus is used 
to constrain the parameters of the plasma model, and the 
background is taken as a region located far from any clus- 
ter emission. We then fit a mekal and internally absorbed 
power-law model to the data (Model II). Note that, other 
more complex absorbers, such as partial covering absorption 
models (zpcfabs, pwab) do not provide good fits. 

MACS J1931.8-2634 - The background-subtracted nu- 
clear spectrum is clearly dominated by non-thermal emis- 
sion; an absorbed (Galactic) mekal model provides a very 
poor fit even if we limit the fitting range to the 0.5 — 2 keV 
energy band. We therefore concentrate on fitting different 
non-thermal models to the spectrum. A simple absorbed 
(Galactic) power-law does not provide a good fit. However, 
a power-law with additional internal absorption at the red- 
shift of the source is able to reproduce the observed spec- 
trum (see Table 131 Model I). We also try fitting a double 
power-law model to the data, but the model is highly un- 
constrained, and more complex absorbers (zpcfabs, pwab) 
are not able to reproduce the s pectrum. The lumin osity we 
find agree well in the results of jEhlert et al.ll201ll ). 

MACS J0947.2-f 7623 - Fitting an absorbed (Galactic -|- 
internal zphabs) power-law to the background-subtracted 
nuclear spectrum does not provide a good fit at ener- 
gies below 1 keV. A different absorbing component at the 
redshift of the source is needed to explain the emission. 
ICavagnolo et all (|201ll ) modelled the nuclear emission us- 
ing a partially covered absorber (zpcfabs), as well as a 
power-law and 2 gaussian lines to account for features seen 
around 1.8 and 3 keV. However, by using their model, we 
were not able to constrain the parameters, especially those 
for the additional gaussian lines. Their model essentially re- 
quires several parameters to be frozen in order to constrain 



the fit. ICavagnolo et al.l (|201ll ) also fitted a distribution of 
partially covered absorbers with the pwab model, but this 
model is also unable to converge to a solution. We therefore 
choose to calculate the flux of the non-thermal component 
in this source by fitting an absorbed (Galactic -|- zpcfabs) 
power-law model to the emission (Model I) . We also include 
a model with an additional gaussian line around 3 keV, but 
keep the width frozen at 100 eV to help constrain the fit 
(Model II). 

MACS J2046. 0-3430 - Fitting an absorbed (Galactic) 
MEKAL model to the background-subtracted spectrum does 
not provide a good fit, and underestimates the emission be- 
yond 2 keV. This suggest that there is a significant non- 
thermal contribution to the spectrum. If we fit a simple 
power-law model to the spectrum with Galactic absorption, 
the resulting index is unnaturally large (F > 4). Including 
internal absorption at the redshift of the source does not im- 
prove the fit, and the internal absorption converges towards 
a null value. Instead, we consider the possibility that there 
is a contribution of both thermal and non-thermal emission 
to the nuclear spectrum, and apply the same method as for 
RXC J1558.3— 1410 to help constrain the fit. In this case, 
the surrounding 2" — 3" annulus is used to constrain the 
thermal plasma parameters within the 1" circular region, 
and the background is taken as a region located far from 
any cluster emission. We consider power-law models with- 
out (Model I) and with (Model II) internal absorption at the 
redshift of the source, but keep the index frozen (F — 1.9) 
to help constrain the fit. 

MACS J0913. 7-1-4056 - The background-subtracted nu- 
clear spectrum, where the background is taken as the sur- 
rounding 2" — 3" annulus, shows three distinct components: 
one at soft X-rays (1 keV) and two at hard X-rays (a bump 
around 3 keV and an emission line at 4.6 keV with a ~60 eV 
width). The location of the emission line coincides with the 
redshifted 6.4 keV Fe K emission complex. A simple ab- 
sorbed (Galactic) power-law with an additional emission line 
does not provide a good fit at soft X-rays, even if we add 
an internal absorption at the redshift the source. However, 
the emission seen around 1 keV can be modelled by adding 
a complex absorber at the redshift of the source (zpcfabs). 
The fit is constrained if we keep both the energy (4.6 keV) 
and width (60 eV) of the line frozen. The resulting values are 
shown in Table [4] (Model I). We also consider the possibility 
that the emission seen around 1 keV is of thermal origin. 
In this case, we model the background-subtracted nuclear 
spectrum as an absorbed (Galactic) mekal and power-law 
with emission line model, where we keep the emission line 
energy and width frozen. We also keep the abundance of 
the thermal component frozen at the best-fitting value for 
the surrounding 2" — 3" annulus. The fit is constrained, but 
the best-fitting power-law index is negative. If we add an 
internal absorption at the redshift of the source (zphabs), 
the model is not able to fit properly the emission seen at 
soft X-rays, and the upper limit of the temperature is not 
constrained. Freezing the temperature at a similar value de- 
rived from a mekal model fitted to the surrounding 2" — 3" 
annulus does not provide a good fit. It is therefore more 
likely that the emission seen around 1 keV is of non-thermal 
origin and we model the 0.5 — 7 keV emission as a double 
power-law, where one of the power-laws is also affected by 
internal absorption at the redshift of the source (Model II). 
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To help constrain the fit, we keep the energy and width of 
the emission line frozen, and first fit the 0.5 — 2 keV energy 
range to constrain the power-law index of the first power- 
law (at soft X-rays). This yields a value of 1.6 ± 0.5. We 
then fit the complex model to the entire 0.5 — 7 keV range 
and keep the index of the first power-law frozen at 1.6. The 
nuclear spectru m of this source has been looked at in de- 
tail by authors (jlwasawa et al.l l200ll : IPiconcelli et al.l l2007l : 
IVignah et al.ll201lh . 

MACS J1411. 3-1-5212 - The nuclear spectrum of this 
source is similar to the one in MACS J0913. 7-1-4056. There 
are two distinct components: one at soft X-rays (~ 1 keV) 
and one at hard X-rays (« 3 keV). A simple absorbed power- 
law cannot account for both components and a more com- 
plex model is needed to explain the emission. We begin by 
fitting an absorbed (Galactic) power-law with a partial cov- 
ering absorber at the redshift of the source (zpcfabs) and 
show the results in Table U (Model I) . The spectrum can 
also be fitted with a double power-law model, where one of 
the power-laws is affected by additional internal absorption 
at the redshift of the source (Model 11). However, we also 
consider the possibility that the emission seen at soft X- 
rays is of thermal emission. We keep the abundance frozen 
at the best-fitting value of a mekal model applied to the 
surrounding 2" — 3" annulus (0.627© ) and let the temper- 
ature, normalization parameters, internal absorbing column 
density and power-law index free to vary (Model III). 

MACS J1423.8-I-2404 - We begin by fitting a mekal 
model to the background-subtracted nuclear spectrum. Al- 
though the fit converges, it slightly underestimates the emis- 
sion seen at hard Xrays (> 3 keV). Instead, we try fitting an 
absorbed (Galactic) power-law to the spectrum (Model I), 
which provides a better fit at hard X-rays. We also consider 
the possibility that there is both a thermal and non-thermal 
contribution to the spectrum, and apply the same method 
used for RXC J1558. 3-1410 to help constrain the fit. This 
method consists of using the surrounding 2" — 3" annulus 
as a proxy for the thermal plasma parameters within the 1" 
circular region. Applying this method, we fit a mekal an 
internally absorbed power-law model to the nuclear spec- 
trum. The background is taken as a background far from 
any cluster emission. For the fit to converge, we also require 
the power-law index to be frozen. In this case, the inter- 
nal absorption converges towards a null value. We therefore 
choose to fit the nuclear spectrum with a simple mekal and 
absorbed power-law model (Model II). 

Cygnus A - The Chandra nuclear spectrum of this 
source shows clearly the presence of a non-thermal compo- 
nent at hard X-rays which can be modelled as an absorbed 
(intrinsic) power-law with an additional gaussian line. How- 
ever, there is also the presence of a faint soft X-ray compo- 
nent that we first model as a power-law. Here, we keep the 
gaussian line energy and width fozen at a value of 6.1 keV 
and of width 0.05 keV to help constrain the fit (Model I). 
We also consider the possibility that the emission seen is 
of thermal nature, and apply the same method as for RXC 
J1558.3— 1410 to help constrain the fit. In this case, the gaus- 
sian line energy and width are also kept fozen at a value of 
6.1 keV and of width 0.05 keV to help constrain the fit 
(Model II). 

Note that in HL2011, we performed a spectral analysis 
for many of the sources with no detectable X-ray nucleus 



such as Abell 1835 and ZwCl 3146 in an attempt to fur- 
ther constrain the upper limit of the non-thermal fiuxes. 
Although we were able to obtain a rough estimate of the 
non-thermal fluxes for these objects, by applying the same 
method as for RXC J1558.3— 1410 (i.e. the method where we 
use the surrounding annulus to estimate the mekal plasma 
parameters of the inner 1" circular region), we stressed that 
most of the parameters in the analysis had to be frozen in 
order to constrain the flt. It was also not clear if the addi- 
tion of the power-law component for these objects actually 
improved the fit, and we proposed that these apparent de- 
tections could be simply due to the manifestation of cooler 
thermal components which is often seen in the central re- 
gions of cool-core clusters. We therefore assume that it is 
reasonable to treat all the derived luminosities for the 19 
objects in our sample with no detectable X-ray nucleus as 
upper limits. 
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